






tlRoimi^^Bflci :' #» iiini^MqP 




1. Report No. 2. Government Accession No. 

NASA TMX -3317 

3, Recipient's Catalog No. 

4. Title and Subtitle 

COMPUTER METHOD FOR DESIGN OF ACOUSTIC 
LINERS FOR TURBOFAN ENGINES 

5. Report Date 

October 1976 

6. Performing Organization Code 

7. Author(s) 

Gene L. Minner and Edward J. Rice 

8. Performing Organization Report No. 

E-8428 

10. Work Unit No. 

505-03 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 


A computer-oriented design procedure for the specification of acoustic liners for turbofans is 
presented. The method begins with an estimate of the noise generated by a turbofan engine. 
This estimate is based upon existing information in the literature. The generated source 
spectrum combined with the results of psychoacoustics leads to a target attenuation spectrum. 
The liner is designed to meet the target attenuation spectrum in a least squares sense by using 
the results from an empirically modified theoretical acoustic propagation analysis, together 
with an empirical impedance model. An example is carried through the discussion in order to 
demonstrate the method. A limited amount of data are presented to illustrate the comparison 


between estimate and experiment. 


17. Key Words {Suggested by Author(s) } 

Acoustic liners; Duct acoustics; Turbofan 
engines; Noise reduction; Perceived noise; 
Liner design; Aircraft noise; Engine noise; 
Internal noise; Liner; Suppressor; Acoustics 


18. Distribution Statement 

Unclassified - unlimited 
STAR Category 71 (rev. ) 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) ^ 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

9.2 

$4. 75 


* For sale by the National Technical Information Service, Springfield. Virginia 22151 























COMPUTER METHOD FOR DESIGN OF ACOUSTIC LINERS 
FOR TURBOFAN ENGINES 
by Gene L Minner and Edward J. Rice 
Lewis Research Center 

SUMMARY 

A design package is presented for the specification of acoustic liners for turbofans. 
An estimate of the noise generation is made based on modifications of existing noise 
correlations, for which the inputs are basic fan aerodynamic design variables. The 
method does not predict multiple pure tones. A target attenuation spectrum is calculated 
that is the difference between the estimated generation spectrum and a flat annoyance- 
weighted goal attenuated spectrum. The target spectrum is combined with a knowledge 
of acoustic liner performance as a function of the liner design variables to specify the 
acoustic design. The liner design method at present is limited to annular duct config- 
urations. The detailed structure of the liner is specified by combining the required im- 
pedance (which is a result of the previous step) with a mathematical model relating im- 
pedance to the detailed structure. The design procedure is developed for a liner con- 
structed of perforated sheet placed over honeycomb backing cavities. 

A sample calculation is carried through the report in order to demonstrate the de- 
sign procedure, and experimental results presented show good agreement with the cal- 
culated results of the method. 


INTRODUCTION 

The reduction of aircraft engine noise by installing acoustically treated surfaces 
inside the engine nacelle has become accepted practice with the advent of current gen- 
eration aircraft. For example there is acoustic treatment in the engine inlet ducts of 
the wide-body transports in production. Such treatment is installed primarily for the 
reduction of tone and broad band noise generated by the fan. As noise regulations for 
new certifications become increasii^ly restrictive, additional suppression will be re- 
quired. One of the ways to accomplish additional suppression is to introduce acoustic- 
ally treated rings in inlets and exhausts. 



This report presents a procedure, in the form of a computer program, to design 
acoustic treatment installations for annular ducts. This design procedure is particularly 
convenient because it permits the methodical examination of alternative choices of 
treatment for noise reduction. Such choices as number of splitter rings traded off 
against lei^h are important because they affect engine nacelle weight, aerodynamic 
losses, and ultimately engine economics. The procedure is one which is based upon 
currently available theoretical and experimental information and thus has relevance at 
the present time but which must be updated and modified as new information and new 
situations arise. It is therefore not claimed to be a final procedure applicable for all 
time. 

The design procedure described combines results from theoretical models of wave 
attenuation in acoustically treated passages with correlations of experimental data from 
full-scale fan suppressors. While theoretical models are satisfactory for laboratory 
duct experiments where the input wave structure is reasonably known and where simple 
geometries are employed, these models have not shown as good an agreement as de- 
sired in full-scale fan and engine experiments. For this reason experimental data from 
fan tests were correla,ted for setting peak attenuation levels and attenuation bandwidths. 
Parametric dependences derived from theoretical results were retained however. 

The theoretical models that exist do not lend themselves directly to the design of a 
liner. The reason for this problem lies in the fact that the theory lends itself to eval- 
uating the acoustic attenuation for a given installation. The general design problem, of 
course, asks for an installation specification to achieve a desired or target attenuation. 
To implement application of the theory to solve the design problem, the theoretical 
analysis was exercised parametrically to cover a wide range of installations. These 
parametric exercises combined with experimental data provide the basis for perform- 
ance maps and equations which are used in the design effort. 

The resulting procedure lends itself very nicely to screening a number of design 
configurations with respect to aerodynamic losses, installation weight, etc. , and thus 
to the selection of the most pragmatic configuration. By introducing experimental data 
correlations from fan and engine tests, the treatment specification program accounts in 
an approximate way for "real- world” effects such as wall boundary layers, duct termi- 
nations, sound modal structure, and other similar effects to the extent that the influence 
of these effects can be generalized for a variety of fans and engines. 


PRCXIRAM APPROACH 

In the section PROCEDURE in this report, methods are presented for estimating 
source noise and for designing a liner to suppress the noise. Such an application is 
shown in figure 1. The source is considered to be a low tip speed fan with the fan jet as 
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the noise floor. The suppressor design methodology employs rectangular duct theory 
and therefore applies specifically to aimular passages such as the ringed inlet or the 
fan exhaust duct of figure 1. Some further work should enable the use of the procedure 
for designing cylindrical- duct suppressors. The suppressor design method could be 
applied to other noise sources having other noise floors. The various parts of the de- 
sign problem fall naturally into a set of distinct, but interrelated, subpackages. As 
such, it is a straightforward matter to implement these subpackages in the form of 
subroutings to be called, as needed by a computer design program. Appendixes A and 
B contain listings of the required computer programs. As new information becomes 
available on the various components of the program, it is possible to incorporate 
changes into the proper subroutine or subroutines. It is also possible to perform the 
necessary steps of the design procedure by hand calculations, which are however 
tedious. 

To design a noise suppressor to meet a required goal, it is necessary to have an 
estimate of the unsuppressed noise source spectrum. There are potentially several 
ways to secure such an estimate, for example through testing of the unsuppressed en- 
gine, or by theoretical prediction. The present procedure employs empirical noise 
correlations taken from the literature and modified somewhat as deemed necessary. It 
is not claimed that this is the best or the only way to achieve the required noise esti- 
mate; rather, it serves as one possible means of getting the needed information. As 
new and better information and/or prediction methods become available, they should be 
substituted for the present methodology. 

A schematic representation of the interrelated parts of the present method is shown 
in figure 2. The design problem begins with an estimate of the noise produced by the 
fan, including the fan exhaust jet noise, which may constitute a floor for available noise 
reduction. This estimate is made using modifications of well-known noise correlations 
from the literature, which require specifications of certain design and operating param- 
eters as inputs. From the calculated noise production and the desired noise goal (i, e, , 
perceived noise level PNL), the required attenuation spectrum can be established. 
Combining the target attenuation spectrum with functional representations of liner be- 
havior in an iterative fashion results in specification of the gross installation charac- 
teristics such as design frequencies of peak attenuation and required acoustic impedance. 
The number of rings, (e.g. , as in fig. 1) and associated passage heights are chosen ex- 
ternal to the computer program. Without the geometric constraints of the installation, 
the goal of these choices would be to obtain the most economically efficient attenuation 
of the noise. 

The required acoustic impedance is used in an approximate mathematical impedance 
model that relates the required impedance to the detailed structure of the liner, includ- 
ing the cavity backing depth, face plate open area ratio, thickness, and perforation size. 
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This step, then, completes the design procedure as it presently exists. 

Figure 2 also indicates the possibility for inclusion of a pressure loss estimate and 
installation weight estimate. Such capabilities would aid in the further optimization of 
liner design, but are not included in the present procedure. 

If a measured fan noise spectrum is available, such a spectrum can be used at the 
appropriate step in the procedure, instead of calculating the fan noise estimate. A 
measured spectrum would be the best obtainable starting point for that specific applica- 
tion. Such a spectrum should be representative of the noise source in flight. An in- 
flight test is one way to measure the source noise. Alternately, if ground tests are 
used, care must be taken to avoid or correct for additional noise generation due to 
ground- static effects which would not be present in flight, such as ground vortices, flow 
distortions, inlet turbulence, and so forth. This point is discussed in further detail in 
reference 1. 

In the course of the development of the liner design procedure, many simplif 3 dng 
assumptions were necessary. This need resulted in part from the very large matrix of 
possible applications and in part from the shortage of available information on certain 
aspects of the problem. For example, expressions for liner bandwidth cannot be ob- 
tained in general for all possible wall construction types because they depend on the 
particular impedance frequency characteristics involved. 

Some of the major specializations and simplifications embodied herein are as fol- 
lows: 

(1) In the noise estimate procedure no multiple pure tones are considered. 

(2) No noise floor other than the fan jet noise floor is considered. 

(3) The acoustic absorber is assumed to be placed in an annular duct with flow 
Mach numbers equal to or less than 0. 6. This restriction arises from the 
theoretical analysis of the duct sound propagation and from the empirical base 
for attenuation rates and wall impedances. It is reasonable to believe that ex- 
tension to other cases could be accomplished if the effects were properly ac- 
counted for, 

(4) The liner itself is of the perforated plate over honeycomb type. Generalization 
to other types of wall construction would require alterations of the bandwidth 
curves and adaptation of the impedance model used for specifyii^ the wall con- 
struction. None of the simplifications contained here appear to invalidate the 
procedure. As additional information is made available, the procedure can be 
changed in a modular sense within the context of the present format. 

With the foregoing discussion providing an overview, the detailed procedure will be 
discussed sequentially in the following sections. 
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ESTABLISHMENT OF TARGET ATTENUATION SPECTRUM 


The noise generation calculation method, as discussed in this section, is presented 
as one means of providing the unsuppressed noise spectra needed for the design of a 
noise suppressor. These procedures are straightforward in their application. How- 
ever, they may not be the best or most direct for all design problems. If the designer 
has information at his disposal which he considers more complete, or superior to that 
presented here, then of course he should use his own information as the preferred 
starting point. In particular, e^erimental measurement of the generated noise spec- 
trum is a good way to secure an unsuppressed noise spectrum. Furthermore, if a de- 
signer has information about noise floors which might limit the noise reduction of a 
suppressor in addition to the fan jet floor discussed here, then that information should 
be taken into account for his purpose. An important point to consider is that the present 
design method provides a versatile framework within which such Improvements can 
easily be made. 

Figure 3 shows an example of a spectrum which has been calculated, for a low 
pressure ratio fan, using the noise correlations discussed in this section. The figure 
is presented here to provide a point of reference for discussion of the different parts of 
the noise estimate which follows. As seen in figure 3, the fan sound spectrum is broken 
into three parts; (1) internally generated tones, (2) internally generated fan broadband, 
and (3) externally generated fan jet noise. The manner in which the composite spectrum 
is assembled will be described in the discussion to follow. 


Prediction of Internally Generated Noise 


The internally generated noise of a low tip speed turbofan has generally been con- 
sidered to be composed of discrete tones associated with the blade passing frequency of 
the fan and of broadband noise related to the turbulent flow of air over the internal sur- 
faces of the machine. From a large number of tests of various engines, these noise 
components were correlated as functions of the fan operating parameters and geometry 
by Smith and House (ref. 2). Another, more recent, noise correlation procedure which 
could be used is given by Heidmann (ref. 3). The main motivation for selection of the 
Smith and House method was its availability at the time the present procedure was or- 
ganized. Either method should provide reasonable noise estimates, and the computer 
could be provided with alternate subroutines to use either or both correlations. 

The Smith and House correlation yields estimates of the noise spectra at the sideline 
angles of maximum noise in both the inlet and exhaust hemispheres. The correlation is 
expressed simply in terms of normally available fan parameters. Comparison of the 
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results of the Smith and House correlation with experimental results of reference 4 and 
with more recent, as yet unpublished, turbofan results has indicated the need to modify 
the Smith and House correlation in order to give reasonably accurate correlations of the 
noted recent results. The Smith and House correlation was developed from machines 
having inlet guide vanes, while the present applications have none. It is assumed that 
the present rotor-stator interaction noise behaves in the same way as the inlet guide 
vane-rotor noise correlated by Smith and House. 

Broadband noise. - For fans with no inlet guide vanes, the correlation for predict- 
ing the sound pressure level at the peak of the broadband noise curve SPLpgg is 

SPLnRR = ’75* + 10. log J + 50. log(-^| + AF (1) 

\ 1000 / 

where Jl is the air flow rate through the fan in pounds per second, Vtip is the tip ro- 
tative speed of the rotor in feet per second, AF is an adjustment for the front- rear 
directivity split of sound from the fan and is a function of the tip relative inlet Mach 
number of the fan. Then, 

AF = 10. log (1 + M^.p), M^.p £ 0. 8 (2) 


In the original method positive and negative values of M^p were used for rear and 
front maximum angles, respectively. However, it was found for the experimental re- 
sults in the Quiet Fan Program that the prediction works well at the front angles but not 
at rear angles where the broadband is about the same as in the front. Therefore aF 
was calculated using positive for both front and rear estimates. Thus the change 

of sign in the original method is ignored for the broadband noise. Note that the calcula- 
tion gives the peak sound pressure level at the front maximum (50°) and rear maximum 
(120°) noise angles, on a 30. 5- meter (100-ft) sideline. The frequency of the broadband 
noise peak in Hertz as given by Smith and House is 


, _ 10 000 
pBB ch 


(3) 


where Ch is taken, for fans without IGV's, to be the average of the axial chords at the 
tips of the rotor and stator given in inches. 

The shape of the broaeband spectrum is estimated by a log normal curve fit of the 
spectral data of SPL - SPLpgg presented by Smith and House (ref. 1) to be 





( 4 ) 
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where 



(5) 


The constant 0. 883 causes a minor adjustment of the actual peak of the bell-shaped 
curve fit of elation (4) from that calculated by equation (3). This shift is not a sig- 
nificant effect, but is required to cause the bell to fit the Smith and House results in a 
best sense. The resultant broadband spectrum is shown in figure 4. 

Discrete tone noise. - The discrete tones are predicted using the following equation: 


SPLx = 85. + 10, 


where Jl is the air flow rate in pounds per second; S/Ch is the ratio of axial rotor- 
stator spacing to mean of the axial projection of the tip chord of the rotor and stator, 

fH 

and n indicates the n harmonic of the blade passing tone. In the case of the tones, 
AF is used as intended by Smith and House in the front and in the rear; that is, both 
negative and positive signs are used on as appropriate. 

The frequencies at which the tones occur are found as follows: 


f = Bn 
60 


( 7 ) 


where RPM is the fan speed in revolutions per minute, B is the number of rotor 
blades, and n is the n*^ harmonic of the blade passing tone (n = 1 for the fundamental). 
Each tone is placed in the one-third octave band containing the frequency of the tone. If 
two or more harmonics fall within a particular band, the intensities are added and the 
resultant sum is used. No attempt is made to divide the power between two adjacent 
one-third octave bands when a tone falls near the theoretical boundary between two 
bands. 

Combined spectra . - Combining the broadband and tone noise leads to the estimated 
spectrum of internally generated noise available for suppression. Recall that no pre- 
dictions of multiple pure tone noise are included in the method; therefore, such cases 
are excluded from consideration. As discussed previously, figure 3 shows an example 
spectrum of this internally generated fan noise combined with the estimated fan exhaust 
jet noise. The fan parameters for the NASA Lewis QF-9 fan (low tip speed and pres- 
sure ratio, ref. 5) were used for the example calculation, and their values are given in 
table I. To assemble the spectrum, the fan broadband noise is used as the starting 
spectrum. At those one-third octave bands where the estimated tone level exceeds the 
broadband, the tone level replaces the broadband point. H the tone level is equal to or 
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TABLE I. - FAN DESIGN PARAMETERS FOR CALCULATION 


OF SOURCE NOISE ESTIMATE OF QF-9 FAN 


Fan speed, RPM 2227 

Number of rotor blades, B .............. 15 

Fan weight flow, Jl, kg/sec (Ib/sec) 404 (890) 

Spacing to chord ratio, S/C. 2 

Tip Mach number, 0.63 

Mean axial chord, Ch, m (in. ) . . . . . . 0.316(12.45) 

Jet exit area, A, (ft^) . 2. 02 (21.76) 

Jet characteristic length (passage 


height, Lj,, m(in.) , . . , , . . . . , , 0. 473 (18, 6) 
Jet velocity, Vj, m/sec (ft/sec) . 177 (579) 


below the broadband, the one^thir'd octave band point remains unchanged. 

There have been some unsatisfactory results with attempts to compare this tone 
noise prediction to measured data for a wide range of fans. Poor comparisons are at- 
tributed in part to the rotor-stator spacing term in the prediction, which may not be 
applicable to fan noise data in static tests. This observation is made because of the 
recently advanced postulate that inlet distortion and inlet turbulence interacting with the 
rotor are major sources of tone noise in static testing of fan stages (refs, 1, 6, 7, 8, 9, 
and 10), If such were the case, it is clear that rotor wake-stator interaction would not 
necessarily be the dominant source of tone noise and that therefore attempts to account 
for the rotor-stator spacing effect could introduce scatter into the comparison between 
predicted and measured tone noise. It is probable that, for fans in which rotor wake- 
stator interaction is the dominant source of tone noise, the rotor-stator spacing de- 
pendence used here is a good approximation. This may apply to the higher harmonics of 
the fan tone statically and to the fundamental tone in flight, provided it is above cutoff. 
For purposes of estimating the spectra of preliminary fan designs for use in examing 
liner tradeoffs, the relations presented are very useful. 


Prediction of Fan Jet Noise 

The externally generated fan exhaust jet noise is important because it represents 
an ultimate noise floor below which the addition of more acoustic treatment would jdeld 
no observable reduction of the far field fan noise and because it contributes to the per- 
ceived noise level. Thus the fan jet noise presents a practical and ultimate limit for 
noise attenuation by acoustic liners, unless the observer might be shielded from the jet 
noise, such as in an engine-over-the-wing configuration. In addition to the fan jet 
noise, there are other sources in a complete turbofan engine which need to be considered 
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as potential floors for such a machine. For example the core engine jet noise, the com- 
bustor noise and the core engine turbine noise may be significant sources which should 
be considered when warranted. Other noise floors, such as self- generated noise by the 
flow over the treatment, may lie above the fan jet noise and thus become the real limit 
in some cases. However, little is known about treatment of self noise. 

In the interest of simplicity, the present development is limited to the specific 
application of a turbofan with its associated single stream exhaust jet noise as the noise 
floor. Extension of the present design method to cover other situations would not re- 
quire significant effort, given the needed information about other noise sources. 

The jet noise is estimated using a procedure based on the "SAE correlation” 

(ref. 11) which has been modified to account for the results reported in reference 12. A 
recently published jet noise prediction method is that of reference 13, which could be 
used as an alternate to that used in this report. The present methodology uses the 
former, since the development of the present procedures preceded the publication of 
reference 13. The level of the jet noise is set by considering the maximum jet noise 
overall sound pressure level OASPL^ at an angle of maximum jet noise on a 30. 5- 
meter (100-ft) sideline as a function of jet size and velocity. Figure 5 shows a plot of 
OASPLj^ - 10 log A as a function of jet velocity. This curve was taken from ref- 
erence 11 and adjusted in level in order to remove jet density from the correlation pa- 
rameter such that the curve passes through the fan data of reference 12. For use in the 
computer, the curve shown in f%ure 5 is represented as 

OASPLj^ = 93. log V. + 10. log A - 176. 2 (8) 

The result of this equation, using v^ in feet per second and A in square feet, is the 
jet OASPL^ at a 30. 5-meter (100-ft) sideline and maximum angle. The jet noise max- 
imum angle in a sideline representation is assumed to occur at 135° from the engine in- 
let centerline. 

In order to estimate the level of the jet noise at the same far field point as the in- 
ternal noise levels (namely 50° and 120° on a 30, 5- meter sideline), adjustments are 
made for the distance and directivity effects. The directivity information is taken from 
reference 12 and is shown in figure 6 as the variation of OASPL directivity index as a 
function of angle. The value of the OASPL at an angle other than the maximum side- 
line jet noise angle can be calculated using the directivity index. The curve of figure 6 
is used to find the difference between the directivity indices at the desired angle and the 
angle of maximum OASPL. This difference, when added to the OASPL^, yields the 
value of OASPL at the new angle on an are passing through the maximum jet noise lo- 
cation (135°, 30. 5-meter sideline). The jet noise estimate is then adjusted to the 50° 
and 120°, 30. 5-meter sideline positions of maximum internal noise by making a distance 
correction, assuming a point source and including a small adjustment for atmospheric 
absorption. 
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The directivity index curve was drawn through the data of reference 12. A table of 
values was taken from this curve and stored in the computer for adjustments to the de- 
sired locations. Since the OASPL directivity and a fixed-shape jet noise spectrum are 
used, it follows that no account is taken of possible variations in directivity of different 
frequency bands. 

The Strouhal spectrum of fan jet noise was taken from reference 11 and was fitted 
with a segmented analytical expression for ease in computation. The curve is shown in 
figure 7. To find the jet noise spectrum for a specific fan, the following procedure was 
used. A dimensionless Strouhal number normalize the spectrum 

where f is frequency, L is the characteristic length, and v. is jet velocity. For 
the present, L was set equal to nozzle annulus height, which was found to work well 
for the spectra of reference 12. To predict one-third octave spectra, a simple level 
adjustment from the octave-band spectrum Is made. 

It is observed that the jet noise correlation, thus modified, may provide a reason- 
able estimate for the fan jet noise from turbofans which are similar in geometry to that 
of reference 12, that is, small annulus height as compared to the diameter and long 
trailing center body. For significantly different fans the results may be in error. This 
comment applies mainly to the use of annulus height in the Strouhal number for place- 
ment of the spectrum, with respect to frequency. 

The total noise picture is formed by energy adding the jet noise spectrum to the in- 
ternally generated spectrum as shown in figure 3. 


Comparison with Measured Noise Data 


Figure 8 shows a sample Comparison at 120° between the data measured on the 
QF-9 fan of reference 5 and the noise estimation results discussed previously. Agree- 
ment is seen to be quite good. The agreement is seen to be poorest at frequencies be- 
low 250 hertz. The discrepancy there is probably due mainly to the fact that there is a 
ground interference problem in the data which is not included in the estimate. The ap- 
parent over-prediction of the blade passing frequency tone (first spike) is due in part to 
the splitting of the measured tone between the 500- and 630-hertz bands. As noted 
earlier such splitting does not occur in the estimate. As a rough comparison, all of the 
measured tone can be placed in the 500-hertz band by adding about 2. 5 decibels there. 
The resulting comparison is reasonably good. It should be noted that the QF-9 results 
were not used in the data base which determined the estimation procedure. 

In contrast with the good agreement shown here, there have been other cases where 
the agreement was not so good. This result is attributed to the general simplicity of 
the correlation procedure in representing a complicated process. It is, however, ar- 
gued that, as long as the fan for which a noise prediction is sought is reasonably similar 
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to those which formed the empirical basis for the estimation procedure, the resulting 
noise estimate should be reasonably good for the purposes of the method. 

Note that no attempt is made to give a noise estimation procedure for fans having 
buzz-saw or multiple pure tone noise. This additional source should be considered 
when appropriate. At present no general multiple pure tone prediction method exists; 
however, an interim method published in reference 3 can be used to calculate an esti- 
mate pure tone noise for high (supersonic) tip speed fans. 


Target Attenuation Spectrum 

The object of placing an acoustic liner in a fan duct is to reduce the annoyance 
(perceived noise) of the noise to an acceptable level. To ascertain the amount of re- 
duction required, the designer needs to have a quantitative measure of the annoyance 
caused by a sound pressure level spectrum such as that of figure 3. Reference 14 pre- 
sents a method relating a quantity called Noys, which is directly proportional to human 
annoyance, to sound pressure level as a function of one- third octave frequency band. 

The effect of this quantity is to provide a weighting function which accounts for the fact 
that some sound frequencies produce greater annoyance than others. Figure 9 shows 
the spectrum of annoyance (Noys) derived from the estimated sound pressure level spec- 
trum of figure 3, using the method of reference 14. It is notable that the frequency 
range from 1000 to 4000 hertz is heavily weighted, and the tone at 500 hertz contributes 
a substantial amount of annoyance. 

Once the target perceived noise level PNL is known, it is possible to find a sound 
pressure level spectrum which meets the goal. A simple first choice for a spectrum 
having the target PNL value is that for which the sound in each one-third octave band 
is equally annoying. Figure 9 includes a flat Noy ’’spectrum” which would lead to 
90 PNdB. This value of Noy was derived using the method of reference 14. The equa- 
tion for perceived noise level is 

PNL = 40 + 33. 21 log n (9) 

where 

i7=0.85Noyjjj^+0.15 £ Noy. (10) 

all 1 /3 oct 

For the case where Noy^ is a constant Noy^, this equation becomes 

77 = (0. 85 + 0. 15 Nj) Noy^. (11) 
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where is the number of frequency bands and Noy^. is the constant value of Noy. 
Rearranging equations (9), (10), and (11) reveals that 


Noy L___ 

^ 0.85 + 0,15Nj 

and 

(PNL-40.)/33.21 

Solving for Noy gives 

(PNL-40.)/33. 21 

Noy = 

® 0.85 + 0.15Nf 


( 12 ) 


(13) 


(14) 


Finally, with this value of Noy, it is possible to enter a table or an equation which 
relates SPL, Noy, and frequency to yield a sound pressure level spectrum with the 
target PNL. A subroutine which employs a curve fit of the relation from reference 14 
is used for the calculation of the spectrum. Figure 10 presents a sound pressure level 
spectrum calculated by this subroutine yielding a constant 7. 2 Noy level and having a 
resulting value of 90 PNdB. 

It should be obvious at this point that there are an infinite number of possible 
sound pressure level spectra which would result in a given goal PNL. The flat Noy 
spectrum is a simple choice for this purpose. The individual designer could, of course, 
make some other choice which could be implemented in a fairly simple way. 

A perceived noise level goal is generally specified at a distance greater than the 
30. 5-meter (100-ft) sideline discussed in the preceding noise estimation procedure. 
Therefore before the difference between the goal spectrum and the noise- generation 
spectrum can be calculated, the generated spectrum must be adjusted to the desired 
distance. In the present case, the generated spectrum is adjusted to the distance of the 
goal, accounting for noise reduction due to distance and excess atmospheric attenuation. 

Once the target spectrum has been determined, it is an easy matter to subtract the 
target spectrum from the predicted noise generation spectrum to obtain the resultant 
target attenuation spectrum for the liner, keeping in mind that the noise cannot be re- 
duced by the liner below the level of the fan jet noise. Figure 11 presents an example of 
a target attenuation spectrum at 120°, using the noise generation estimate shown in fig- 
ure 3 for the 1. 2 pressure ratio (QF-9) fan at design condition and the target flat Noy 
spectrum of figure 10 calculated to reach the 90 PNdB goal. 

In order to establish a target attenuation spectrum, it was necessary to gamine the 
spectrum of sound pressure level at the angle of maximum annoyance and to compare 
this spectrum with a predetermined desired attenuated spectrum. In the previous illus- 
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trative example, this angle of maximum annoyance was 120°, Recall also that the for- 
ward angle of maximum annoyance using the noise estimation procedure of this report 
occurs at 50°. For the situation where the designer haf other information at Ms dis- 
posal on maximum annoyance spectra, he should, of course, use that information. 

The determination of noise estimates at the forward maximum and rear maximum 
angles, aloi^ with noise goals, as discussed preidously, allows the calculation of 
sound attenuation requirements for inlet and exhaust suppressors, respectively. The 
target sound pressure level attenuation spectra are assumed to convert directly to tar- 
get power level attenuation spectra for the suppressors. This assumption implies that 
the inlet suppressor attenuates uniformly over the inlet hemisphere and that the exhaust 
suppressor attenuates uniformly over the exhaust hemisphere. Furthermore, it is 
implied that there is no significant noise effect in the inlet hemisphere from exhaust 
noise, and vice versa. These assumptions are necessary in the absence of sufficient 
information on attenuation directivity effects. It is known, however, that tMs procedure 
should jdeld attenuation estimates for the liner which are conservative at the angle of 
maximum noise. 

Appendix A contains a listing and a discussion of a FORTRAN IV computer program 
to perform the calculations indicated for determination of target attenuation spectra. 

The example problem carried through tMs report is also included as input-output 
listings. 


SPECIFICATION OF LINER DESIGN 

The procedures which are set forth herein apply specifically to the design of an 
acoustic liner for an annular passage having opposing walls treated with identical ma- 
terials. The liner itself, seen schematically in the insert in figure 1, has a perforated 
face sheet backed by honeycomb with an impervious wall behind the honeycomb. Any 
deviations from this type of liner would have to be taken into account by altering some 
details of the design procedure presented here. 

In effect, the present procedure considers the design of a sir^le annular acoust- 
ical passage, and it must be assumed that any parallel passages (for example when 
splitters are present, as seen in the inlet of a research fan, fig. 1) will be constructed 
to accomplish the same attenuations. Alternatively the passage considered most im- 
portant for design might be the outermost passage. The profile of sound intensity, as 
well as area distribution, often is Mgh at the outer wall and falls off toward the center, 
thus concentrating a large portion of the sound energy in the outer passage. Neverthe- 
less, the inner passages caimot be ignored in the overall suppressor design, because 
they will allow a noise leak if they are poorly designed. Furthermore, the passage 
heights and number of splitters must be compatible with the overall geometry of the 
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installation. The number of splitters and related passage heights are not presently 
handled automatically by the program, but are set externally in an iterative manner in 
such a way as to accomplish the cjesired attenuation within the physical constraints of 
the duct, while allowing for splitter thickness. 

The liner design proceeds in essentially two phases. The first phase employs the 
combined e3q>erimental results and theoretical prop^ation analysis in order to specify 
the required acoustic design; that is liner peak attenuation frequency, liner length, and 
optimum liner impedance. In the second phase of the design, the tuning and optimum 
impedance requirements of the first phase are translated into hardware designs using 
an empirical mathematical model for wall impedance to specify perforated sheet open 
area ratio and honeycomb backing depth. Such items as facing sheet thickness, hole 
diameter, and honeycomb cell size are selected as inputs to the program based on 
availability of materials, ease of fa.brication, cost, or other pragmatic considerations. 
These parameters do have an influence on acoustic performance and cannot be changed 
arbitrarily once the remainder of the installation has been specified based on some 
choice of these parameters. 


Phase 1 - Acoustic Design 


The parameters to be specified in this part of the design procedure are passage 
heights, tuning frequencies, length of treated surface and optimum impedance at the 
frequency of peak attenuation. A number of reports in the literature, such as refer- 
ence 15, 16, and 17, have presented results of analytical calculations of liner attenu- 
ations. Some typical kinds of problems considered have been propagation of plane wave, 
or spinning modes, or least attenuated modes. For purposes of the present design 
procedure in order to limit the complexity of parametric relations, the curves of at- 
tenuation for least attenuated mode in a rectangular duct were selected as presented in 
reference 17. The rectangular duct is a good approximation to the annular passages in 
an inlet with splitter rings and in the fan exhaust. The items to be considered in this 
section involve the prediction of the liner performance in terms of the variables dis- 
cussed in the preceding paragraph. The basis for these predictions are theoretical 
analyses of acoustic propagation for both optimum and off-optimum impedances. To 
improve agreement with experimental observations, an empirical constant has been 
introduced to adjust the levels of the estimates. Establishing the spectral estimates by 
this means as functions of the noted parameters enables the design choice of the param- 
eters to achieve a target attenuation. 

Peak attenuation at optimum impedance . - One of the results of the propagation 
analysis is the prediction of sound attenuation at the frequency of peak suppression at 
optimum impedance, for the rectangular duct least- attenuated mode. This peak sup- 
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pression was predicted as a function of duct Mach number M, duct length L, duct 
height H, peak frequency fp, and speed of sound c. For the least attenuated mode 
these parameters were combined; and the prediction was presented concisely in terms 
of dimensionless parameters as shown in figure 12. 

Experimental data at the Lewis Research Center have shown that the attenuation at- 
tained in practice for ringed liners is less than that predicted by the theoretical curves 
of reference 17. Therefore the curves presented here in figure 12 have been adjusted 
to show suppression levels only 0.4 times those in reference 17, 

It is normal practice to summon a large collection of data in order to specify an 
empirical factor. In the present case experimental results from tests on a single inlet 
liner were used to arrive at the factor of 0. 4. This was done because of the extreme 
difficulty of attempting to correlate measured attenuations from suppressors which were 
off optimum by varying amounts. It is, however, a relatively easy exercise to employ 
the prediction procedure, once the empirical factor has been determined, to calculate 
estimates for other suppressors and thereby verify the empiricism against other data. 
Limited comparisons of this nature are presented later in this report. The agreement 
has been good for cases examined to date. However, such comparisons should continue 
in order to further build confidence in the empiricism. 

The details of the technique used to specify the 0.4 constant factor were as follows. 
Spectral attenuation data from reference 18 for the inlet suppressor in the QF-IA fan 
testing program were used. The details of the acoustic design for this liner are given 
in table II. 


TABLE n. - QF-IA LINER USED FOR EMPIRICAL 

ADJUSTMENT OF PEAK ATTENUATION CURVES 

[Duct Mach number, M, -0. 32; duct temperature, T, 
288 K (519° R); static pressure, p, 9. 32x10^ N/m^ 
(13. 5 Ib/ln. ^); duct height, H, 20 cm (7.9 In. ); sheet 
thickness, t, 0. 051 cm (0. 020 in. ); hole diameter, d, 
0. 081 cm (0. 032 in. ).] 


Ogen 

Back 

Equivalent 

Calculated 

Calculated 

area 

cavity 

length, 

peak 

resistance 

ratio. 

depth, 

L 

frequency 

ratio, 

percent 

1 




V 

®/«opt 

cm 

in. 

cm 

in. 

Hz 


2.5 

0. 51 

0, 20 

38 

15. 

2741 

1.4 

2,5 

2,23 

, 88 

48 

19. 

1184 

4.6 
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The attenuation spectrum, rather than the attenuation at a single frequency, was used in 
order to increase the data base for the factor. Making the calculated estimate for this 
spectrum, required considerations of off-optimum impedance and liner bandwidth be- 
havior which to this point have not been discussed here. The use of these details was 
consistent with the discussions to follow. 

Figure 13 presents the comparison between the QF-IA data and the adjusted pre- 
dicted estimate. Note that the fit is quite good. Included among the experimental re- 
sults are attenuation spectra from several different speeds. These data with the ex- 
ception of the tones, exhibit a good correlation about a single faired curve. 

The reason for the large shift between the pure theoretical prediction and the ex- 
perimental results (as reflected in the factor 0. 4) is not fully understood. A possible 
cause may lie in boundary layer and mode shape effects, which are not explicitly in- 
cluded in the present method and which could therefore tend to de -optimize the imped- 
ance, thus reducing attenuations below those e 3 q)ected. 

As seen in figure 13, there are several data points (tailed) which lie above the gen- 
eral trend. These points represent attenuation of tones which may damp rapidly be- 
cause they occur as spinning modes, as discussed in reference 16 or as higher order 
radial modes. The present procedure does not consider spiiming modes or higher order 
radial modes, and therefore these data points were not included in setting the value of 
the empirical constant. If tones damp faster than broadband sound as these data sug- 
gest, a parallel method for predicting the attenuation of tones is needed. However, the 
present procedure when applied to tones would provide conservative estimates of atten- 
uation. 

For use in the computer, the curves of figure 12 were e 3 q>ressed in the form of 
equations. The entire curve for zero Mach number was fitted with a single equation for 
which the form was motivated by the shape of the curve as follows; 



where B’s are curve-fit constants, given in table Ill(a), and rj is the dimensionless 
frequency parameter 77 = Hfp/c. 

The curves at M / 0 are fit by considering the deviation from the M = 0 curve, 
such that 


L 

H 




+ 0. 4 Delj^ 


( 16 ) 


16 


M=0 



where 


D, 


Delj^=D^, 


k2 


- D, 


>k4 + ^ 


k3 


(17) 


The k subscript indicates that the curve-fit constant D’s are different at different 
Mach numbers. The values erf are listed in table ni(b) for the Mach numbers 
shown in figure 12. For other Mach numbers the peak attenuation is found by interpol- 
ating. This is done in the computer by using a pair of interdependent interpolation 
routines, taken from references 19 and 20. At values of Mach number outside the range 
shown in table HI, the use of extrapolation may not be reliable. 


TABLE in. - CONSTANTS FOR CURVE FIT OF 
LINER PEAK ATTENUATION (EQ. 15) 

(a) No flow case (dact Mach number = 0) 


=0.41792 
Bg = 0. 57188 
Bg = 0. 64086 
B^ = 0. 23905 
Bg = 14. 167 


(b) Duct flow case (Mach number general), eq. (17) 


Mach 

number, 

M 

Integer 

k 

Integer i 

1 ^ 

1 

2 

3 

4 

Values ^ 

-0.4 

1 

1.5475 

45. 497 

18. 811 

2.4012 

- .2 

2 

1.1661 

10. 937 

9. 3986 

1. 1438 

0. 

3 

0. 

0. 

0. 

0. 

.25 

4 

-2. 0987 

30.741 

12.255 

2.4699 

.4 

5 

-3. 5123 

49. 095 

13. 861 

3.5251 


6 

-4. 1881 

48. 918 

14. 121 

3.4267 


There are instances when it may be advantageous to design a liner with an impedance 
which is off-optimum at the frequency of peak attenuation. Such a design could have 
the advantage of a broader bandwidth than an optimized liner. However, when the liner 
resistance is not optimized, the magnitude of peak attenuation is adversely affected. 

Performance at off-optimum impedance. - This section presents background in- 
formation which led to the design procedure presented in other sections. The reading of 
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this section is not necessary to the application of the method. In order to estimate an 
attenuation spectrum for the liner to be designed, it is necessary to know how the liner 
behaves at frequencies away from the frequency of peak attenuation. Therefore the 
bandwidth characteristics of the liner must be determined. The peak sound power at- 
tenuations shown in figure 12 and the liner optimum impedance study reported in ref- 
erence 16 provide information at only one frequency, the peak attenuation frequency, 
and this only if the liner is optimized. For the rest of the attenuation spectrum, off- 
optimum performance must be considered. 

Determination of the attenuation spectrum is made by combining results from the 
theoretical propagation analysis with calculations of impedance from a semi-empirical 
impedance model for a specific kind of wall material. The propagation analysis re- 
sults in sets of contours of constant attenuation on an impedance plane plot. These con- 
tours are generated under some idealized assumptions. It is assumed that, at any given 
frequency, the impedance, specific resistance 0, and specific reactance x (which are 
parameters normalized by pc, the free space impedance) can be varied arbitrarily. 
Thereby a relation between resistance and reactance, which leads to constant values of 
attenuation, can be found. Changing the frequency results in a new set of contours. The 
impedance of an actual liner has a specific dependence on frequency, on liner geometry, 
and on environmental conditions. This dependence is approximated by the impedance 
model. The combination of the propagation analysis and the impedance model to give an 
attenuation spectrum is a somewhat complex process which is described in greater de- 
tail in reference 21. A brief illustration of the basic ideas involved will be presented 
here. 

In reference 21 it was found that the constant attenuation contours could be approx- 
imated by a series of circles in the impedance plane described previously. This series 
of circles was represented by a single equation containing only the maximum possible 
sound power attenuation and optimum wall impedance as inputs. Figure 14 shows an 
illustrative example using conceptual contours at three different frequencies: £^, 0. 5, 
and 2 f. (solid, dashed, and dash-dot curves, respectively). On any one of the contours 
the attenuation is a constant for that frequency. The attenuation shown on the curves are 
all relative to the peak attenuation, which is taken as unity, at f.. The peak attenuations 
at 0, 5 f . and 2 f. were set relative to that at f. by using the inverse liner dependence 
of attenuation on the frequency parameter as shown in figure 12 for values of 77 > 2. At 
one-half the intermediate frequency, the peak attenuation is two; while at twice the in- 
termediate frequency, the peak attenuation is one-half. The optimum resistance 
and reactance Xqp^ vary as 

^ CX £ (18) 

^opt J 
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which is valid for M = 0. If M ^ 0, for equation (18) to be valid, the mode must be 
well cut-on (propagating), as was shown in reference 16. 

As noted previously, the complex value of impedance was allowed to vary arbitrar- 
ily in order to construct the contours of figure 14, The impedance of an actual liner 
does not vary in an arbitrary way. In order to use the contours of figure 14 to calculate 
attenuations, it is necessary to introduce the actual behavior as a function of frequency. 
To demonstrate the technique, two separate illustrative examples will be considered. 

In the. following discussion of impedance, it is convenient to consider the resistance 
and reactive components individually. The actual resistance will be assumed independ- 
ent of frequency; this is a reasonable assumption, in particular whenever there is sig- 
nificant flow past the liner surface. The actual wall specific reactance is modeled by 
the following cotangent relation: 


^ = -coti!S. =-cotlj!a (19) 

c c 

where co is radian frequency, b is cavity backing depth, f is frequency (Hz), and c 
is the speed of sound. For the present illustrative examples, equation (19) will be ap- 
proximated by 

X oc - i (20) 

The first example is a liner with its values of resistance and reactance at frequency 
f . coinciding with the optimum values for that frequency; This point is shown as an 
open diamond in figure 14. At other frequencies the liner resistance is a constant, as 
assumed previously, and the reactance is described by equation (20), Therefore the 
liner impedance value lies on a horizontal line passing through the optimum point. Two 
other points on this line are called out: the open square at a frequency of 0. 5 f^ and the 
open triangle at a frequency of 2 fj. When 0 = 1 and x = "1 wall impedance coin- 
cides with optimum impedance for f^, and the relative attenuation is one. ForO. S f^^ 
the actual wall specific resistance is d = 1, and the specific reactance is x = and 
this point must be associated with the dashed contours. The relative attenuation is thus 
seen to be 9. 53. For 2 fj the actual wall specific resistance is 6 - 1, and the specific 
reactance is x = ~0. 5, which is associated with the dash-dot contours. The relative 
attenuation is seen to be slightly less than 0. 2, say 0. 19. Therefore the fall off in at- 
tenuation with frequency is seen to be faster on the high frequency side than on the low 
frequency side of the optimum. This is seen to be due to the displacement of the actual 
wall impedance and of the optimum loci with frequency, as well as to the reduction of 
the maximum possible attenuation with increasii^ frequency. 

The second example is a liner with a value of resistance which is too large to give 
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maximum possible attenuation at any frequency, as shown in figure 14 as solid symbols. 
The liner characteristics are chosen such that at frequency f j, the actual liner reactance 
equals the optimum value of reactance for that frequency (diamond- shaped symbol). At 
other frequencies the liner resistance and reactance are specified as before. The values 
at 0. 5 fj and 2 f^ are illustrated by the square and triangles, respectively. For the 
overdamped (large resistance d) off-optimum liner the peak attenuation is seen to be 
around 0. 56 (solid curves). The one-half and double frequency points are seen to have 
relative attenuations of about 0. 51 and 0. 34, respectively. Although the peak attenuation 
is down to one-half the optimum peak, the attenuation spectrum is seen to be relatively 
flat, which is typical of overdamped systems. The two examples discussed here have 
shown in principle the way that attenuation is calculated for a liner with known or as- 
sumed impedance behavior. By considerations of this nature the complete liner behav- 
ior as a function of frequency (referred to as bandwidth) can be calculated for a liner 
with optimum or off-optimum values of impedance. The general relations to follow de- 
pend on the precise nature of the actual liner impedance, and any changes from the as- 
sumptions used here would necessarily alter the general relations. 

Off-optimum peak attenuation. - For a liner which is not optimized, the frequency 
of peak attenuation occurs at the frequency where the reactance equals the optimum re- 
actance. The peak value of attenuation in the spectrum is reduced in amplitude from that 
possible with an optimized liner because the resistance does not equal the required op- 
timum resistance at this frequency. 

The complete determination of liner off-optimum performance required considering 
a large number of cases similar to those in the preceding section; however, the more 
precise formulation of equation (19) was used instead of the approximation of equa- 
tion ( 20 ) to develop the relations to follow. This section presents the results for off- 
optimum peak attenuation, and the next section presents the results for bandwidth, which 
also arises from off-optimum impedance values. 

The calculated ratio of peak attenuation to maximum possible attenuation is shown 
in figure 15 as a function of liner resistance ratio. This quantity is also a function of 
the frequency parameter 77 (for over damped liners); thus there is a separate curve for 
every value of 77 . As e^qpected, the maximum peak attenuation is obtained when the 
liner is optimized (0/0 Qp|. = 1 ), either overdamping (0/0Qpt > 1 ) underdamping 
(0 /0Qp|. < 1) reduces the peak attenuation. As will be seen in the following sections, 
one might design an overdamped system to increase the attenuation bandwidth at the ex- 
pense of the peak attenuation. 

Liner bandwidth . - The previous illustrative example, in connection with figure 14, 
demonstrated the principle of the method used to calculate the spectrum of attenuation 
for a liner. The procedure discussed could be used for every value of frequency of in- 
terest. This procedure has been simplified such that the attenuation spectrum is as- 
sumed to be characterized by the sketch in figure 16. Once the peak attenuation is es- 
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tablished, the only remaining information necessary for an estimate of the attenuation 
spectrum are the half attenuation frequencies above and below the peak. The frequen- 
cies fj and fg are the frequencies at which one-half of the peak attenuation is ob- 
tained. This representation of the attenuation spectrum should be adequate for prelim- 
inary design calculations. 

In figure 17 the bandwidth frequencies fj and f^ are shown as functions of design 
parameter 77 for a liner with an optimized resistance. The values of f^ and fg have 
been calculated from a combination of the equal attenuation contours and the liner actual 
impedance behavior. The bandwidth is seen to be extremely small for low rj with some 
broadening at high 77. This is due to the stretching of the equal attenuation contours of 
figure 14 at higher values of the frequency parameter (resistance also higher). It should 
be recalled from the illustrative example that, as the bandwidth increases, the peak at- 
tenuation decreases with increasing frequency parameter (see fig. 12). Therefore a 
careful examination of both phenomena must be made in order to design a liner. 

The low-frequency half attenuation point fj/fp for off- optimum liners is shown in 
figure 18(a) for several values of design frequency parameter. For underdamped liners 
( 0/0 opt < l)j ^i/^p essentially constant. Thus very little gain in bandwidth occurs; 
but pe^ attenuation is lost, as can be seen in figure 15. For an overdamped liner 
( 0 / 0 opt > 1) considerable bandwidth increases are obtained with increasing resistance. 

The high-frequency half attenuation point ^ 2 ^^^ for off-optimum liners is shown in 
figure 18(b) with 77 as a parameter. Note that, as with fj/fp, f 2 /^p remains constant 
for ^/^opt ^^’ and as increases above one, the bandwidth increases. The 

curves of figure 18(b) at values of 77 < 10 are adjusted from those one would derive 
directly from use of the cotangent relation of equation (19). In brief the adjustment was 
required to eliminate unrealistic behavior of the cotai^ent function at certain values of 
its argument, namely 0, tt, and so forth. This point is discussed in greater detail in 
reference 21. 

For purposes of the present design procedure, the curves of figures 15 to 18 are 
used to calculate bandwidth and off- optimum liner performance estimates. The curves 
in these figures were fit by a series of straight line segments on the log-log basis for 
ease of implementation in the computer. For values of v which fall between these 
curves, interpolation subroutines from references 19 and 20 were used to calculate the 
needed values. 


Application to Liner Design 

Procedure. - Once the peak attenuation and the bandwidth have been specified as in- 
dicated, the matter of designing the liner is accomplished by choosing design variables 
so that the goal attenuation spectrum is fit in a least-squares sense. This procedure can 
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be summarized functionally as follows: The total estimated liner attenuation spectrum 
is the sum of the spectra from the N liner segments given by 

AdB(fc) = £ AdB.(fp., L^, H, y (21) 

j=l 

where f„ is one-third octave center frequency, j is an index for a particular liner 
segment, fpj is the peak frequency for that liner segment, L. is the lei^h, and H is 
the duct height. The sum of the squares of the deviations between the estimated liner 
attenuation and the target attenuation spectrum over the one-third octave center fre- 
quencies of interest is 

Y = ^ [AdB(y - AdB,(y]^ (22) 

The extent of the one-third octave spectrum used in equation (22) is somewhat arbitrary. 
Selection of this range must be based on the shape of the target spectrum. 

The least-squares calculation makes the selection of fp’s and L’s in such a way as 
to minimize Y. The number of peak frequencies N and the duct height H are set by 
practical considerations of the geometry. The minimization of Y is done by selecting 
each variable so that each of the derivatives 

u = — - 0 (23) 

9w 

where w represents each one of the variables fp. and Lj. Thus the deviation is 
minimized with respect to each variable. The derivatives are evaluated by a numerical 
procedure in which each of the design parameters is guessed and then incrementally 
stepped. New values of parameters are chosen by Newton-Raphson iteration 

in which 


w 

new 


= w 


old 



(24) 


The iteration stops when w^^^^ is within an acceptable tolerance of for each of 

the variables sought. 

Example results . - Usii^ the techniques discussed previously, an example liner was 
designed for the target attenuation spectrum of figure 11 for the QF-9 fan aft duct. The 
resultii^ liner specifications are shown in table IV together with some input environmen- 
tal parameters listed in the headnote. These parameters were assumed rather than 
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actual for this machine. This example is also included in appendix B to demonstrate 
the application of the liner design program. 


TABLE IV. - SAMPLE ACOUSTIC LINER DESIGN PARAMETERS FOR QF-9 


EXHAUST DUCT SUPPRESSOR 


[Duct Mach number, M, 0. 4^ duct temperature, T, 288 K (519® R); static pressure, 
p, 1. 034x10® N/m (15 psia); duct height, H, 0. 246 m (9.70 in.).J 


Section 

Peak 

Dimensionless 

Length to 

specific 

Specific 


frequency, fp, Hz 

frequency, 

height ratio, 

resistance, ^ 

reactance, ^ 


V 

L/H 

e 

X 

1 

502. 

0.36 


*^0. 595 

■1 

2 

- ■ - ■ ■ 

3011. 

2. 18 


1.115 



^Normalized by pc (free space impedance). 
^Off optimum by choice 6=2 


The liner design given in table IV has two sections with frequencies of peak atten- 
uation of 502 and 3011 hertz. The number was chosen based on the shape of the target 
spectrum. The program calculated the particular frequencies. The value of H used 
here was chosen in order to achieve low rj and large L/H and thus increased liner 
effectiveness. This duct height could be achieved by placing a single splitter ring in the 
existing QF-9 duct and expanding the outer casing by the thickness of the splitter, 

A design estimate attenuation spectrum as well as the target attenuation for this 
liner is presented in figure 19. At the lower frequency tuning, a somewhat broader 
banded liner than the optimal case was considered advantageous in the interest of broad- 
ening the bandwidth without prohibitive length penalties. Therefore the resistance ratio 
at the peak frequency was selected at 0/®Qpt; = 2. This, however, had an adverse ef- 
fect on the peak attenuation. The magnitudes of the peak and of the bandwidth were de- 
termined from figures 15 to 18. 

Recall that in the design procedure it is assumed that the surfaces which face each 
other in the duct are treated with identical materials. There can, however, arise 
practical considerations, both from a flow and a layout point of view which favor the 
staggering of materials. For example, as in the present design illustration, a low fre- 
quency tuning requires a large cavity depth, which would lead to a large splitter thick- 
ness. Thick splitters decrease inlet pressure recovery, and they are thus undesirable. 
Therefore it is most convenient to place such low frequency treatment on the outside and 
inside walls, and not on the splitter. Calculations for such a case have shown that as 
long as opposing walls are treated at frequencies which are reasonably close, there will 
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be no significant differences from the ease with opposing walls the same. Proper ac- 
count must be taken of the required total area of each acoustic liner, as calculated by 
the analysis. Placing treatment on only one wall of a passage, in effect, would require 
approximately doubling the length of that treatment, for example. One is, in effect, 
moving the treatment from the splitter to the outer and inner walls. The placement of 
low frequency treatment on only outer and iimer walls can be done only if there is either 
one or no splitters because if there were more than one splitter, there would occur a 
pass^e without any low frequency treatment. 


Phase 2 - Mechanical Design 

The value of qptimum impedance required to give the predicted sound attenuation is 
calculated using the theoretical propagation model of Rice (ref, 15). Once the required 
value of acoustic liner impedance (e.g., optimum or off -optimum per choice) has been 
specified, the mathematical model relating impedance to the geometric parameters can 
be used to calculate open area ratio cr and backing depth b. The basic equations of the 
impedance model which are semiempirical were reported by Groeneweg in reference 22. 
However, these equations did not include the effect of the free stream steady flow on the 
impedance. In the present model empirical terms are added to account for this effect 
on specific resistance and reactance. The details of this mean flow effect are taken 
from reference 23, The following resistance equation is used to calculate the open area 
ratio a: 

j3 +l/^^ + 4eM 

a«_L li! (25) 

20 


where 




(26) 

(27) 

(28) 


and where t is facing plate thickness, d is hole diameter in the facing plate, v is 
kinemamatic viscosity, w is circular frequency at the liner peak attenuation frequency 
(o) = 27rfp), p is air density, is flowby Mach number (negative in the inlet and 
positive in the exhaust duct), and p^ ^ is the standard reference pressure used in cal- 
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culating sound pressure level (2xl0"^ Mbar), and OASPL is the duct overall sound 
pressure level. A reasonable estimate for the OASPL in-duct can be computed by 
assuming that acoustic energy is conserved between the duct and the far field. The far 
field energy can be assumed distributed uniformly over the duct cross section and the 
OASPL can be computed thereby. 

The cavity depth (backing depth) b is found by solving the following reactance eq- 
uation: 


b = £ arctan 22_ (29) 

0) (t + 6)a> - x<7C 

where 6 is an inertial end correction for the air in the hole of the facing sheet; that is, 

;,_ 0. 85d (l. - 0.7V^) (3QJ 

1. + 305. |Mp 

where is the duct free stream Mach number. 

In the previous equations most of the terms are evaluated by a knowledge of flow 
conditions in the duct. Quantities which as yet appear unspecified are sheet thickness 
t and hole diameter d; there are no ’’equations*' to calculate them. They are there- 
fore set by alternate considerations, Such as the nature of the environment, tuning fre- 
quency, availability of materials, and ease of fabrication. The nature of installation or 
environment (e. g. , high temperature flow) may set a minimum thickness t. The tuning 
frequency, if very low, may require a very thick facing plate or even inset tubes, as 
used in reference 24, in order to keep the backing depth reasonably small so the sup- 
pressor will fit into available space and/or minimize aerodynamic pressure losses if 
used on a splitter ring. 

Table V presents the open areas and backing depths, as well as the input sheet 
thicknesses and hole sizes for a liner in continuation of the example of the previous 
sections and described in table IV. Liner section 1 has a large backing depth, hole di- 


TABLE V. - LINER MECHANICAL DESIGN, EXAMPLE LINER 
OF TABLE IV, FOR QF-9 DUCT 


Section 

Sheet 

Hole 

Open area 

Back cavity 


thickness. 

diameter^ 

ratio, 

depth. 


t 


c 

1 

(T, percent 

1 

) 


cm 

in. 

cm 

in. 


cm 

in. 

1 

3.8 

1.5 

0.64 

0.25 

23.7 

7.06 

2.779 

2 

.051 

. 020 

.127 

.050 

12.0 

1.37 

0.541 
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ameter, and sheet thickness because of the low frequency of peak attenuation and the 
fact that exhaust ducts (M > 0) generally require larger backing depths than inlets re- 
quire- The large thickness can be simulated with reasonably light we^ht welding 
tubes of the required length to the back side of the facii^ sheet as was done in refer- 
ence 24. This problem does not occur in liner section 2- 

Appendix B contains a listing and a discussion of the liner design FORTRAN IV 
program. Also tabulated with the listii^ is a sample output of the program for the ^- 
ample problem carried throi^h the report. The t 3 rpical run time required for the pro- 
gram is less than a minute on the IBM 7094 computer. 


Design Chart Based on Propagation and Impedance Models 

The optimum impedance computed by the propagation model can be combined with 
the impedance model to give design charts as a function of fp, tj, environmental param- 
eters such as duct static pressure p and temperature T, and other geometric vari- 
ables of the liner. This design chart is not intended to replace the computations re- 
quired by the preceding design methodology. Rather, the chart is an illustrative ex- 
ample of the trends contained within the calculation procedure. Of course, to cover all 
cases of interest, many charts would be needed. 

An example of such charts is presented in f^re 20 where the back-cavity depth is 
shown as a function of dimensionless frequency parameter rj at selected values of fre- 
quency fp, and the open area ratio a is also given. Each of the parameters tabulated 
there also has a determinant effect on the chart, thus a new chart would be required if 
any one of these parameters were changed. As a point of reference, the locus of points 
for a constant H of 0. 305 meter (1 ft) is included in figure 20(a). This locus is de- 
termined by coimecting the points on the various curves for b at different frequencies, 
notii^ the interdependence of H, fp, and t} through the definition of 77 . All points to 
the right of the constant H curve have values of H > 0. 305 meter (1 ft), and all points 
to the left have values of H < 0. 305 meter (1 ft). This curve emphasizes the point that 
fp and 7 j are not independent if H is specified. 

The curves of figure 20(b) show that, over a wide frequency range, required open 
area ratio does not vary significantly, once frequency has been accounted for throv^h 
the frequency parameter tj. The tolerance on the value of a after fabrication is cer- 
tainly much wider than the variation of design cr with f_ in the curves shown. 

X7 


Comparison of Calculated Liner Performance with Experimental Results 
An idea of how well the present design procedure works in estimatii^ attenuation 
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can be gained by comparing an attenuation estimate with some recent experimental re- 
sults. At present there have not yet been tested any suppressors which were designed 
in accordance with the present procedures. Therefore it is necessary to use the data 
which are available from an existing design and make a performance estimate for the 
case tested. Recall that a set of unpublished data from the QF-IA program was used to 
specify the empirical constant in the design procedure. Obviously these data should be 
in good agreement with estimates made by the design procedure . This point was dem- 
onstrated by figure 13. Other data using the same inlet liner, but mounted on different 
fans, have shown similarly good agreement. As pointed out in an earlier section, fur- 
ther comparisons are needed in order to establish more confidence in the estimating 
procedure used here. The inlet had treated splitter rings in place, similar to those 
shown in figure 1 on another fan. Obviously the engine had no centerline drive shaft. 

In order to make the performance estimate for this inlet, the duct formed by the outer 
wall and the outermost ring was analyzed. Thereby the assumption was made that the 
inner passages were performing at least as well as the outer one and that the behavior 
in the outer passage dominated the far field results. The description of this outer pas- 
sage is given in table VI. As seen in the table, the suppressor had two different sec- 
tions, and each section was calculated to be off-optimum by the present design proce- 
dure. 

Selected for consideration is a set of data (not related to the examples discussed 
elsewhere in this report) from the General Electric Quiet Engine C Program. The en- 
gine and nacelle are described in reference 25; however, the data used here are un- 
published results from tests performed at the NASA Lewis Research Center, 


TABLE VI, - SPECIFICATIONS OF GENERAL 


ELECTRIC, ENGINE C INLET LINER FOR FKHJRE 21 


[Duct Mach number, M, -0. 4; duct temperature, T, 
288 K (519° R); static pressure, p, 9, 66X10^ N/m^ 
(14, Ib/in. ^); duct height, H, 15 cm (5, 9 in,); sheet 
thickness, t, 0, 081 cm (0, 032 in, ); hole diameter, d, 
0,127 cm (0,050 in,)] 


Open 

area 

ratio, 

percent 

Back 

cavity 

depth, 

b 

Length, 

L 

Peak 

frequency, 

V 

Hz 

Resistance 

ratio, 

®/^opt 

cm 

in. 

cm 

in. 

10 

7 

2.54 

.76 

1.0 

.30 

51. 

43. 

20 

17 

1064 

2385 

3.44 

.81 
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The data were taken in such a way as to examine the effect of the inlet liner on at- 
tenuation at takeoff engine speed, forward maximum noise angle, and 61 -meter (200-ft) 
sideline. The experimental attenuation spectrum was found by comparing the SPL spec- 
tra for the hard inlet baseline and the fully treated inlet. Both configurations had a 
full aft duct suppressor such that aft-radiated internal sovind should not be a contributor 
at the forward angle of maximum noise. 

The comparison between the measured attenuation and the estimated performance 
curve is shown in figure 21. The comparison is seen to be excellent, except at very 
low frequencies and at the tone, as was the case in figure 13. 

While the data and the predicted curves presented in figures 13 and 21 show excel- 
lent agreement, further comparisons of this nature are required in order to gain more 
confidence in the prediction method. 


CONCLUDING REMARKS 

This report presents a design package for acoustic liners for turbofan engines. 

The procedure is comprised of an estimate of the noise output of the fan, a goal noise 
spectrum and attendant target attenuation spectrum, and finally a liner specification 
based upon the calculated target spectrum. The procedure is limited to the fan alone 
and does not consider other noise sources ihxthe target attenuation estimate. 

There are no experimental results yet available for a liner designed according to 
the present design procedure. However, a rough comparison was made between a 
performance estimate and measured attenuation for an engine inlet suppressor. In view 
of the fact that the measured data were probably approaching a noise floor, and the lin- 
er was off optimum according to the present model, the comparison is reasonable. 

The design procedure for acoustic liners was developed from results of theoretical 
studies of wave propagation in acoustically treated ducts and from experimental re- 
sults of full-scale fan and engine tests. The procedure applies only for annular flow 
passages as opposed to open cylindrical ducts. This limitation is a result of the fact 
that experimental data on cylindrical ducts are not well correlated by current theoretical 
models. The probable cause is spinning modes which come into play in cylindrical 
ducts and which have been analyzed only recently. Further work is needed to reconcile 
these recent analyses with measured results. Considerable effort would probably be in- 
volved, but the payoff would be an expansion of the design procedure to handle cylindrical 
ducts. 

Because experimental results are incorporated into the procedure, it can be as- 
sumed that a number of real-world effects are included that are neglected or approxi- 
mated in most theoretical studies. Additional comparisons between measured data and 
the correlation used in the present procedure are needed to verify and refine the em- 
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pirical adjustment to the theoretical peak attenuation curves. 

As presented, the procedure has its greatest utility in generating a number of can- 
didate suppressor designs useful for tradeoff evaluations. A selected final design would 
be subjected to considerations that are less generalized than those presented in this 
procedure. It has been pointed out that updated behavior for all of the various factors 
can be incorporated as they are needed or become available. This is an important con- 
sideration because it is recognized that the design of suppressors is relatively in its 
infancy. It is inevitable that new and better information will become available. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, Ai^st 19, 1975, 

505-03. 
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APPENDIX A 


NOISE ESTIMATING PROGRAM DESCRIPTION 
Description 

This computer program estimates the noise generation by a turbofan and uses it 
to establish target attenuation spectra required to reach a set noise target. Input data 
consist of basic geometric and flow information about the fan, as well as atmospheric 
conditions, and the noise target information. The printed and punched output results 
include spectral data for the various components of the noise estimate and the goal at- 
tenuation- spectra at forward and rear angles. In adjusting the separate parts of the 
estimate to a desired location, spherical divergence, and atmospheric absorption are 
taken into account. 

The main program and required subroutines are listed in this appendix, using 
FORTRAN IV computer language. The required data input cards are itemized in de- 
tail as comment cards at the beginning of the main program. Therefore, specific input 
requirements will not be discussed here. 

The printed output resulting from this program is given in tables Vn to XVin. The 
numbers refer to example printout pages included after the main program and sub- 
routines. The example shown is for the case carried throughout the main body of the 
report. 

This program also produces punched output which is compatible with the liner 
specification program input requirements. These punched cards include the variables 
which are explained on comment cards in the liner specification program, appendix B: 
two throwaway explanatory cards, ISTART, ISTOP, GOAL(I), NF, SPL(I), SLJET(I), 
COND(I), and OASPL for both front and rear angles, The cards are arranged with all 
of the previous information punched for the front and then for the rear in block form. 
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Main Program and Subroutines 


$IBT TZ WHLTNG 
C 

:: THIS PROSR&M ESTia&TES FAN NOISE SPECIRA AT 50 AND 120 OESREES AND 

C DErERMINES GOAL, ATTENUATION SPECTRA TO REACH A TARGBI PNL. 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 

c 


BASIC FAN OPERATING PARAHETERS^ DISTANCE, AND THE PNL TARGET 
ABE REQUIRED FOR THESE ESTIMATES. 

A PREDICTED SPECTRUM CONSISTS OF DISCRETE TONES AND BROAD BAND NOISE 
PREDICTED VIA A MODIFIED SMITH AND HOUSE PROCEDURE AND OF JET NOISE PREDI 
CTED VIA A MODIFIED SAE PROCSDJRE. 

DATA REQUIRED 

XPLN1 (I) 12A6 

CARD TO EXPLAIN FAN BEING PREDICTED 
XPLN2(I) 12A6 

SECOND CARD TO EXPLAIN FAN BEING PREDICTED 
FLOW, VR, SOC, RPM, ROTORB, AMC 6F10.0 

FLOW = FAN FLOW RATE, LB/3EC 
VR = FAN-AIR TIP RELATIVE VELOCITY, FT/SEC 

FOR THE CASE OF NO IGV»3, THIS IS ASSUMED TO BE TIP ROTATIVE SPEED. 
SOC = ROTOR-STATOR AXIAL GAP (SPACING) DIVIDED BY ROTOR TIP 

AXIAL CHORD. 

RPM = FAN ROT&TIVE SPEED 
ROTORB = NUMBER OF ROTOR BLADES 

AMC = ROTOR- STATOR MEAN AXIAL TIP CHORD, INCHES 
VJ, AJ, AH 3F10.0 

VJ = JET VELOCITY, FT/SEC 
AJ = NOZZLE AREA, SQUARE FEET 

AH = NOZZLE CHARACTERISTIC LENGTH (ANNULUS HEIGHr), INCHES 
DKIND, AKIND 10X,A1,9X,A5 

DKIND = BLAMX or S, FOR R\DIUS OR SIDELINE CALCULATIONS, RESPECTIVELY 
AKIND = BLANK OR NOATM, FOR ATM ATTN OR NOT, RESPECTIVELY 
DIST, T, RH 3F10.0 

DIST = RADIOS OR SIDELINE DISTANCE, FEET 
T = TEMP, DEGREES F, FOR ATM ATTN. 

RH = RELATIVE HUMIDITY, PERCENT, FOR ATM ATTN. 

AF, AR 2F10.0 

AF = APPROX INLET FLOW AREA, SQUARE FEET. 

AS = APPROX AFT DUCT AREA, SQUARE FEET. 

THESE VARIABLES ARE USED FOR CALCOLATION OF IN-DJCT OA5PLS, AND 
NEED BE SPECIFIED ONLY APPROXIM ATELY FOR SUFFICIENT ACCURACY. 

PNL FI 0.0 

PNL = TARGET PERCEIVED NOISE LEVEL AT CONDITION SPECIFIED BY 
DKIND AND DIST 

THE ABOVE SET OF DATA IS FOR ONE FAN CALCULATION, FOR A NEW FAN A SINGLE 
COMPUTER RUN CAN BE USED, BY ADDING ANOTHER SET OF DATA CARDS. 

SUBROUTINES REQUIRED 
TONES 
BROAD 
JET 
GOAL 
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C ATMSPR 

C BASPAT 

C ATMAT 

r SIDLAT 

C TOGTHH 

C COMPAE 

C PNDBEH 

:: XPLANB 

c 

COMMON /FRBQ/ F(27) 

COMMON /CONSrS/. NP,RPM,S0TORB,FLOW,VH,AMACHF,AHACHR,3OC 

COMMON /COND/ DKIND, AKIND, T, RH, DIS P 

COMMON /PARTS/ STONE (27 , 2) ,SBaOi 0 (27, 2) ,SJET(27,2( 

COMMON /MANY/ N,IT0NE(27) 

COMMON /CHORD/ AMC 

COMMON /DBLF/ AMACHM 

COMMON /MESSG/ AMESSG 

COMMON /ANGLES/ ANG1,ANG2 

COMMON /XPLAN/ XPLN 1 (1 5) , XPLN2 ( 1 5) 

DIMENSION SGOAL (27) ,SATTN (27,2) ,SL(27,2) 

DIMENSION Fa(27) 

DIMENSION ANGL (2 ) 

NAMELIST /NHV FLOW, FR, SOC,RPM, 8 OTORB, AMC ,7J, AJ, AH, DKIND, DIST, AKIN 
1D,T,RH,PNL,AMACHF,AMACHH 

DATA (FD(I) ,Is1, 27) /50. , 63 . , 80. , 1 00. , 1 25. ,160. ,200. ,2 50. , 3 15. , 400. 

1.500.. 630. .800.. 1000. . 1 250 .. 1600. . 2000. .2500. . 3150. . 4000. . 5000. .63 

200. . 8000.. 10000.. 12500. .16000.. 2 0000./ 

DATA AMESS1 ,AMESS2 /6HrOTAL ,6H JET/ 

DATA DTEST, ATEST /1HS,5fl / 

DATA SIDEL,RAD /4HSIDE, 6 HRADIU3 / 

DATA FRNT,8ER /SHFRONT, 4HREAR/ 

NF*24 

DO 1 1=1,27 

1 P(I)=FO(I) 

ANG1=50. 

AN32=120. 

ANGL (1) =ANG1 
ANGL (2) =ANG2 

1035 CONTINOE 

READ (5,2) (XPLNI (1) ,1=1,12) 

2 FORMAT (12A6) 

READ (5,2) (XPLN2 (I) ,1=1, 12) 

READ (5,10) FLOS, VR, SOC, RPM, ROTORB, AMC 

10 FORMAT (6F10.0) 

BEAD (5,11) VJ,AJ,AH 

11 FORMAT (3F10.0) 

READ (5,12) DKrND,AKIND 

12 FORMAT (10X,A1 ,9X,A5) 

READ (5,109) DIST, T,8H 

109 FORMAT (3F10.0) 

READ (5,1090) AF,AR 
1090 FORMAT (2F10.0) 

READ (5,110) PNL 
113 FORMAT (F10.0) 
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n ( j n 


END OF INPOT DATS READING 

WRITE (6,1000) 

10)3 FORtlAT (1H1) 

WRITE (6,120) 

12) FORaAT (60X, 10HINP0T DATA//) 

WRITE (6,1200) (XPLN1 (I) ,1=1, 12) 

WRITE (6,1200) (XPtN2(t) ,1=1,12) 

1230 FORMAT (35X,12a6//) 

WRITE (6,121) FLOW,yR,SOC,RPB,R3TORB,AMC 

121 FORMAT (35X,6F10. 1//) 

WRITE (6,122) VJ,AJ,AH 

122 FORMAT (35X,3F1 0,2//) 

WRITE (6,123) DKIHD,AKIND 

123 FORMAT (45X,A1 ,9X,A5//) 

WRITE (6,1230) DIST,T,HH 

1210 FORMAT (35X,3F1 0. 1 //) 

WRITE (6,1231) AF,AR 

1231 FORMAT (35X,2F1 0. 1 //) 

WRITE (6,124) PNL 

124 FORMAT (35X,F10,2) 

WRITE (6,1000) 

AMArHM=VR/1116. 

C 

~ CALCOIATB INDIVIDOAL SPECTROM COMPONENTS. 

C 

CALL TONES (N, ITONE, STONE) 

CALL BROAD (SBROAD) 

CALL JET (AJ,AH,VJ,SJET) 

C 

C NDW ASSEMBLE THE COMPONENTS AND ADJUST TD DESIRED SIDELINE OR aAOIOS. 

C 

CALL TOGTHR (NF , NA, S L, OAS PLF, OA3 PLR) 

AMESSG=AMESS1 
CALL ATMSPR(SL) 

DO 15 J=1,2 

CALL PNDBER (1 , NF, SL (1 , J) , PNDB1 , PNDB2) 

GO TO (13,14) , J 

13 PNDBF=PMDB1 
GO TO 15 

14 PNDBR=PNDB1 

15 CONTINUE 
WRITE (6,101) 

101 FORMAT (1H0//////////) 

WRITE (6,102) 

102 FORMAT (SOX, 25HPREDTCTSD HARD WALL PNL*S//) 

CALL XPLANE 

WRITE (6,103) PNDBF 

103 FORMAT (60X,8HPNDBF = ,F5. 1//) 

WRITE (6,104) PNDBR 

104 FORMAT (60X,8HPNDBR = ,F5.1//) 

IF (DKIND .ED. DTBST) DOMD=SIDEL 
IP (DKIND .NE, DTE3T) DOMD=RAD 
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HRITE (6,1040) DXST,DtI"ID 
1040 PORaar (55X,3HAT ,F5.0,5H feet, IX, A6) 

WRITE (6,1000) 

&BESSG=ABESS2 
CALL ATMSPfi(SJET) 

CALL GOAL (PHL, NF, SGOAL) 

CALL COMPAS (HP, PBL, DI5T,SGOAL, 3 JET, 3L, SATPH) 

CALCULATE IN-DUCT OASPL'S 

PI=3. 1415927 
8STD=100. 

R1=BSTD/SIN (PI*50. /180. ) 

R2=RSTD/SIN (PI*60./180.) 

205 CONTINUE 

ODF=OASPLr+10.*ALOG10 (4.*PI*B1*ai/AF) 
ODB=OASPLR+10. ♦ALOG10 (4. ♦PI*E2*S2/AS) 

WBITE (6,1001) 

1031 P ORB AT (1H0////) 

WBITE (6,206) ODF,ODR 

205 FORBAT (55X, 17HIN-DUCT OASPLF = , FI 0.1/ 
163X,9HOASPLR = ,F13.1) 

WRITE (6, 1000) 

DO 210 J=1,2 
FRER-FRNT 

IF (J .EQ. 2) FRER=RER 

PUNCH 2060, FREB, (XPLN 1 (I) ,1=1 , 1 2) 

2050 FORBAT (13a6) 

PUNCH 2060, (XPLH2(T) ,1=1,12) 

ISTART=1 
ISrOP=NF 
DO 2063 I=1,NF 
I=I 

IF (I START .GT. 1) GO TO 2061 

IF (SATTN(I,J) -GE. 0. .AND. SATTN(IH,J) . GT. 0.) 
ITEST=ISTAET+3 

2051 CONTINUE 

IF (SATTN(I,J) .GT. 0. .AND. SArrH(I+1,J) .LE. 0.) 
IF (ISTOP .LE. ITEST) ISTOP=HF 

IF (SATTN(I,0) .LT, 0. .AND. I .LT. ISTOP) ISTART 
2053 CONTINUE 

PUNCH 2065, ISTART, ISTOP 
2 065 FORBAT (213 ,54 X, 12HIST ART, IS TOP) 

PUNCH 208, (SATTN(1,J) ,I=1START, ISTOP) 

PUNCH 207, NF 
207 FORBAT (I3,57X, 2HNF) 

PUNCH 208, (SL(I,J) ,I=1,NF) 

203 FORBAT (8F10. 1) 

PUNCH 208, (SJET (I,J) ,I=1,NF) 

PUNCH 209, DIST, DUHD,ANGL{J) 

209 FORBAT (10X,3HAr ,F10.0,5H FEET, IX, A6, IX, F10.0, IX, 
GO TO (2090,2093), J 
2 090 PUNCH 2091 , ODF 
2091 FORBAT (F10. 1, 50X,5HOASPL) 

GO TO 210 

2093 PUNCH 2091, ODR 
213 CONTINUE 
N0CHN3=1 

IF (NOCHNG .EQ. 1) GO TO 1005 

STOP 

END 


ISrAPT=I 

ISTOP=I 
= 1 


7HDE3REES) 
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$IB?r3 TONE 

SOBBOOTINE TONES (H ,ITONE , STONE) 

C CALCOLftTE TONES FOB BPF AND ITS HARMONIOS, ACCOBDINS TO METHOD OF 

0 SMITH AND HOUSE 

C 

C PAR A WE TEES 

C N NUMBER OF ONE THIRD OCTAVE BANDS C0NTAININ3 TONES, 

C ITONE(I) WHICH ONE THIRD OCTAVE BANDS CONTAIN TONES. 

C STONE SOUND PRESSURES IN TONES AT VARIOHS FREaUENCIES. 

COMMON /FREQ/ F{27) 

COMMON /CONSTS/ NF, RPM, HOTORB, PL OW, V R, AM ACHF, AM ACHR, S OC 

COMMON /DELF/ AHACHN 

DIMENSION ITONE (27) , STONE (27,2) 

FF=RPH*ROTORB/60, 

AN=F(NF)/FF 

N=AN+1. 

AN=N 

FTEST=AN*FF 
FACTR=2.**.1667 
FGOOD=FACTR*F (NF) 

IF (FTEST ,GT. F300D) N=N-1 

IHARM=0 

ICNT=0 

TBM1=10.*ALOG10 (FLOW) 

TRM2=50.*ALOG10 (VR/1000.) 

TRMF=10.*ALOG10 ( 1. -AHACHN) 

TRHR=1 0.*ALOG10 (1. 4-AHACHN) 

TRH3=-10.*ALOG10 (SOC) 

TBM3=-20.*ALOG10 (SOC) 

PART=85. +TRH 1 ♦THM2+TRM 3 

DO 30 1=10, NF 

1=1 

FLO=P (1-1) *FACTR 
FUP=F(I) *FACTR 
PWRF=0. 

PHRR=0. 

ID UM=0 

19 AN=IHARM*1 
FTEST=AN*FF 

IF (FTEST .GT. FLO .AND. FTEST . LT. FOP) 30 TO 20 
IF (IDUM) 30,30,21 

20 IHARM=IHARH+1 
IDUM=1DDM+1 
TRM4=-20.*ALOG10 (AN) 

SLF=PARTtTRMF+TRM4 

SLR=PART+TRMR+TRH4 

PHRF = PWRF4- 10.** (SLF/10.) 

PWBR=PWRR+10.*«= (SLR/10.) 

GO TO 19 

21 ICNT=ICNT+1 

STONE(ICNT, 1) =10.*ALOG10(PWRF) 

STONE(ICNT,2) =10. *ALOG10 (PWRR) 

ITONE (ICNT) =I 
30 CONTINUE 
N=ICNT 
C 

C PRINT OUT TONE SPECTRA. 

C 


35 



O U O U f I u 


HBITE (6,101) 

101 POSH AT (1H0) 

WHITE (6,100) 

105 FOH9AT (47X,36HPAH BLADE PASSIH3 TOKE ASD HARHOHICS//) 

WHITE (6,110) 

115 FOBHAT (44X,43HAT 100 FT. SIDELINE, AND 50 AND 120 DEGaBBS//) 
CALL XPLANE 
WRITE (6,102) 

102 FORMAT (53X,30H FREQ STONBF STONER //) 

DO 105 1=1, H 

I=I 

IFHQ=I 

IFB=ITONE(I) 

105 WRITE (6,106) F (IFR) ,STONE (IFR2, 1) ,SrOHE (IFRQ, 2) 
lOi FORMAT (50X,F10.0,2P10. 1) 

WRITE (6,1000) 

1 050 FORMAT (1H1) 

RETURN 

END 


SIBFTC BROADS 

SOBRODTINE BROAD (SBROAD) 

COMMON /FREQ/ F(27) 

COMMON /CONSTS/ NF, RPM, ROTORB, FLOW, 7 R, AMACHF, AMACHR,S5C 
DIMENSION SBROAD (27,2) 

COMMON /CHORD/ AMC 
COMMON /DELF/ AHACHN 
DIMENSION B(4) 

DATA (B(I),I=1,3) /50-514, -88348, .99493/ 

C THIS SOBROOTINE COMPUTES THE BROAD BAND NOISE SPECTROM FOR A SINGLE 

C STAGE FAN ACCORDING TO THE SMITH AND HOUSE PROCEDURE 

C WITH THE EXCEPTION THAT THE DiaECTIONAL PROPERTIES 

C ARB ASSUMED TO BE GIVEN BOTH FRONT AND HEAR BI THE FRONT END FACPOR. 

C 

C ARGUMENT 

C SBROAD THE CALCULATED BROAD BAND SPECTRUM 

FPK= 10000. /AMC 

SPLPK=75.*10.*ALOG10 (FLOW) +55. *ALOG10 (VR/1000.) 

TRMF=10.»ALOG10 (1.-AMACHN) 


TRMF2=10.*ALOG10(1.+AMACHM) 

SPLPK2=SPLPK+TRMF2 


SPLPK=SPLPK+TRMF 

SIGMA=B(3) 

A=B(1) 

FRC=B(2) 
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ra ca o n <a n <3 o o ta n 


ROH GENERAtE SPECTRUII 

DO 30 r=1,NE 

PR=F (I) /FPK 

ARG=ALOG10(FR/FRC) 

AHG=ARG/SIGMA 

AEG=ABG*ARG 

FIT=A*EXP(-ARG) 

DECR=FIT-50. 

SDOilsSPLPK+DECR 

S BROAD (I, 1)-SDUa 


3D0M2=SPLPK2+DECR 
25 SBROAD (1,2) =SD0H 
25 5BR0AD{I,2) =SD0H2 


30 CONTINUE 

PRIST OUT FAN BROAD BAND SPECrRA. 

WRITE (6,101) 

101 FORMAT (1H0) 

WRITE (6,100) 

WRITE (6,110) 

113 FORMAT (!*4X,43HAT 100 FT. SIDELINE, AND 50 AND 120 DEGREES//) 
103 FORMAT (50X,29HFAN BROAD BAND NOISE SPECTRnM//| 

CALL XPLANE 
WRITE (6,102) 

102 FORMAT (58X,20H FREQ SBROAD //) 

DO 105 1=1, SF 

105 WRITE (6,106) F (I) ,S3R0AD (I, 1) 

105 FORMAT (55X, F10.0,F10. 1) 

WRITE (6,1000) 

1033 FORMAT (1H1) 

RETURN 

END 


SIBFTC GOALL 

SUBROUTINE GOAL (PNL, NP, SGOAL) 

C CALCULATES A SPECTRUM, SGOAL, WRICH WILL HELD A CONSTANT SOY CURVE 

C AS A FUNCTION OF FREQUENCY. 

C THIS IS ONE OBVIOUS EASY WAX TO SHOOT FOR A DESIRED PNL. 

C 

C ARGUMENTS 

C PNL = DESIRED TARGET PERCEIVED NOISE LEVEL, PNDB. 

C NF = NUMBER OF FREQUENCY POINTS IN SPECTRUM 

C SGOAL = SPECTRUM GIVING DESIRED TARGET PNL 


37 



COMMON /FREQ/ F(27) 

DIMENSION SGOAL (27) 

ANF=NF 

EX= (PNL-40.) /33. 21 
S0MNOy=10.**EX 
ftNOYC=SUMNOY/(.85 + . 15*ANF) 

DO 30 1=1, NF 

I=I 

AI=I 

IF (I .GT. 9) GO TO 25 

DB=(ALOG(ANOYC) /ALDG{10.) )* (22.«-AI) -2.5«AI*-64. 

TEST=85.-2.*AI 

IF (DB ,LT. TEST) SO TO 30 

DB64=100.+13.-1.4*AI 

GO TO 26 

25 rBH=.04*(AI-19.5) *(AI-19.5) 

DB64 = 100.-11.*EXP {-TRM* (1.+TRM) ) 

26 DB=(10./(ALOG (10.) *ALOG10 (2.) ) ) *AL03 (ANOYC/64.1 »-DB64 
30 SG0AL(I)=DB 

WRITE (6,101) 

101 FORMAT (1H0) 

WRITE (6,100) PHL 

10} FORMAT (45X,31HeONSTANT NOY SPECTROM TO REACH ,F5.1,5H PMDB//) 
WHITE (6,110) ANOYC 

113 FORMAT (55X, 15HC0NSTANT NOY = ,F5.1//) 

WRITE (6,102) 

102 FORMAT (5RX,20H FREQ SGOAL //) 

DO 105 1=1, NF 

105 WRITE (6,106) F ( I) , SGOAL (I) 

106 FORMAT (55X , FI 0. 0 , F 10. 1 ) 

WRITE (6,1000) 

1 003 FORMAT (1 HI) 

RETURN 

END 


$IBFTC ATMS 

SUBROUTINE ATMSPH(SL) 

C 

C CALCULATE NOISE LEFELS AT A DESIRED RADIUS OR SIDELINE, 

C ACCOUNTING FOR TN7ERSE SQUARE AND ATMOSPHERIC ATTENUATION. 

C ASSUMES SL INPUT SPECTRA ARE AT A SIDELINE, WITH ATMOSPHERIC 

C ADJUSTED SPECTRA MAY BE AT DIST SIDELINE OR HADrUS 

C may or may not have atmosphere ACCOUNTED FOR IN ADJUSTING FSOM RSTO 

C TO DIST. 

C OKIND AND AKIND ARE THE DETERMINING CONTROLS. 

C 

C SUBROUTINES REQUIRED 

C SIDLAT 

C BASPAT 

C ATM AT 

COMMON /FREQ/ F(27) 

COMMON /COND/ DKIN D, AKIND ,T, R H, 0 I ST 
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COHMON /CONSTS/ NF, RP!1, H0T08B, FLOH,¥B, AMACHP, AM ACHR, S3C 
COMMON /MESSG/ AMESSG 

DIMENSION AI (2) ,H (2) , SINAI (2) , IF A (27, 2) ,SL(27,2) 

DIMENSION ATA(27) 

DIMENSION BOLD (2) 

DATA DTEST,ATEST /1 HS, 5HNOATM/ 

DATA AI (1),AI(2) /53 . , 1 20. /, NM/2/, NB/3/ 

DATA ANSD, ANSA /6HB ADIBS , 3HATM/ 

DATA SIDBL /4HSIDE/ 

PI=3. 1415927 
BSTD=100. 

SIB50=SIN(PI*50./180.) 

SI H6 0= SI H ( PI *6 0 . /1 3 0 . ) 

B (1) =DIST/SIB50 
H(2)=DIST/S1B60 

IF (DKIND .EQ. DTBST) GO TO 10 
B(1)=DIST 
B(2)=DIST 
10 CONTINOE 

BOLD (1)=BSTD/SIN50 
BOLD (2) =BSTD/STN60 
24 IF (AKIND .EQ. ATEST) GO TO 22 
DO 17 J=1,NM 

17 CALL BASPAT (SL (1 , J) , BSTD, NP, HB, I, BH,IFA { 1, J) ,B ( 3) , APA) 

DEL50=B (1) -RSTD/SXN50 
DBL60=B (2) -BSTD/SIN60 
DO 20 1=1, NF 

SL (I, 1) =SL (I, 1) -TFA (1,1) *DEL50/1 000. 

2 0 SL(I,2) =SL (I ,2) -TFA (1,2) *DEL60/1000. 

22 CONTINOE 

ADJUST LEVELS TO NEW SIDELINE 08 BADIOS, INVEBSE SQOARE. 

DO 27 1=1, NP 
DO 27 J=1,NH 
27 SL(I,J)=SL(I,J) -20.*ALOG10(R(J) /ROLD(J) ) 

26 CONTINOE 

DO 23 1=1, If F 
DO 23 J=1,NM 

IF (SL(I,J) .LT. 0.) SL(I,J)=0. 

23 CONTINOE 

PBINT OOT ATMOSPHEBIC ATTENOAPION SPECTBA. 

WRITE (6,100) 

10) FOBMAT (46X,37HATM3SPHEaiC ATTENOATION PER 1000 FEET) 

BBITE (6, 101) 

101 FORMAT (1H0) 

WRITE (6,110) 

11) FORMAT (44X,45HACCOONTING FOB SPECTBOM SHAPE, FRONT AND BEAR//) 
WRITE (6,102) 

102 FORMAT (53X,30H FREQ DELDBF DELDBR //) 

DO 105 1=1, NF 

1 05 WBITB (6 ,106) F (I) , TFA (1,1) ,T FA (I, 2) 

105 FORMAr (50X,F10.0, 2F10. 1) 

WBITE (6,1000) 

10)) FOBMAT (1H1) 


PRINT OOT PREDICTED SPECTRA 
WRITE (6,201) AMESSG 
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201 FOBMAT (62X,A6//) 

aSItE (6,200) DIST 

203 POBHAT (40X,21 HPBEDICTED SPECfBA AF ,F5.0,24H FT., 50 AHD 120 OSSB 

1EES//) 

IF (OKI HD .BQ. OTEST) OOtID=SIDBL 

IF (DKIN0 .HE. DTEST) DOMD=ANSD 

IF (AKIND .EQ. ATEST) D0nA=ATE3T 

IF (AKIHD .HE. ATEST) DOHA=AHSi 

HRITE (6,202) D0flD,D0f!A 

202 FOBEAT (35X,21HDISTAHCE 30HDITIDH = ,A6,25H, ATHOSPHEBE CONDITION 

1= ,A6//) 

CALL X PLANE 
HBITB (6,203) 

203 FOBMAT (53X,30H FREQ SPLF SPLB //) 

DO 205 1=1 ,NF 

205 HaiTB ( 6 , 206 ) F ( I) , SL (I, 1) ,SL (1 , 2) 

206 FOBMAT (50X,F10. 0,2F 10. 1) 

HBITE (6,1000) 

BE TORN 
END 


$IBFTC T03TH 

SUBBODTIHE TOGTHB ( HF, HA, SL, OASPLF, OASPLR) 
ASSEMBLE THE VARIDHS PARTS OF THE SPECTBA. 


COMMON /PARTS/ STOH E (27 , 2) , SBB 3A D (27 , 2) , S JBT (2 7 , 2) 
COMMON /FBEQ/ F(27) 

COMMON /MANY/ N, ITONE (27) ,NF 
DIMENSION SL(27,2) 

NA=2 

DO 33 J=1,NA 

NCNT=1 

S0M=0. 

DO 30 1=1 ,NF 
1=1 

N»HICH=ITONE (NCNl) 

SND1=SBROAD(I, J) 

IF (I .EQ. NWHICH) GO TO 23 
GO TO 25 

2 3 SNDT=ST0NE(NCNT,.7) 

NCNT=NCNT+1 

IF (SNDT .GT. SND1) 3ND1=SHDr 
2 5 SN02=SJET (I, J) 

C 

2 SOM IS OSED IN OASPL CALCOLATION. 

e 

SUM=SOM4-10. ** (SNDl/13. ) 

ASND=10.** (SHDl/10.) +10.** (SND2/10.) 

SND=10. *ALOG10 (ASND) 

30 SL(I,J)=SND 

GO TO (31,32) , J 

31 OASPLF=1 0.*ALOG10 (SOM) 

GO TO 33 
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32 OASPLB=10.*AiOG10 (SOM) 

33 COHTISaE 

PaiNT OOT COMBIHED SPECTRA HT 100 FEET SEDELISE. 


WRITE (6,101) 

01 FORMAT (1H0) 

WHITE (6,100) 

0) FORHAr (35X,61HC0MBIMED FAN SPECTRA AT 100 FEET SIDELIHE, 50 AND 1 
120 DEGREES//) 

CALL X PLANE 

WHITE (6,102) OASPLF,OASPLR 

02 FORMAT (37X,35HINTERSALLT GENERATED OASPL FRONT = , F5. 1 ,5X, 7 HREAR 
1= ,F5.1//) 

WRITE (6,103) 

03 FORMAT (53X,30H FREQ 5PLF SPLR //) 

DO 105 1=1, NF 

05 WRITE (6,106) F (I) ,SL (I, 1) ,SL (I, 2) 

105 FORHAT (SOX, F10. 0, 2F10. 1) 

WRITE (6,1000) 

033 FORHAT (1H1) 

RETDBN 

END 


BFTC COM PA 

SOBROOTINE COHPAR (NF, PNL,DIST,S30AL,SJET,SL,SATrH) 

COMPARE PREDICTED SPECTRA HIT3 TARGET SPECTROH AND 
DETERMINE REQOIRED ATTENOATION SPBCTRA. 

COMMON /FREQ/ F(27) 

DIMENSION SG0AL(27) ,SL(27,2) ,SJBT(27,2) ,SATTN(27,2» 

DO 25 J=1,2 
DO 25 1=1, NF 

IP (SSOAL(I) .LT. SJBT(I,J)) GO TO 23 
SOFF=SGOAL(I) 

GO TO 25 

23 SOFF=SJBT(I,J) 

2 5 SATTN(I,J)=SL(I,J)-SOFF 
WRITE (6,101) 

101 FORHAT (1H0) 

WRITE (6,100) 

103 FORMAT (45X, 37HATTEN0ATI0N SPECTRA REQOIRED TO REACH//) 
WRITE (6,102) PNL,DIST 

102 FORMAT (48X,F5.1,9H PNDB AT ,P5.0,14H FEET SIDELINE//) 
CALL XPLANE 

WRITE (6,103) 

103 FORMAT (53X,30H FREQ SATTNF SATPNR //) 

DO 105 1=1, NP 

105 WRITE (6,106) F (I) ,SATTN (1,1) ,3 ATTN (1,2) 

105 FORMAT (SOX, F10. 0, 2F10. 1) 

RETORN 

END 
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$IB?TC PNDB8B 

SOBSOOTINE PNDBER (I 1 , 12, DB,PNDB 1 , PNDB2) 

C THIS SUBBOnTINE COBPUTES THE PUL FOB SPECrROH DB, WHICH IS AH ARHAI WITH 

: THE INDEX I CORRESPONDING TO THE NnHBERS OF THE THIRD OCTAVE BANDS. 

C 

C ARGUMENTS 

C II = FIRST BAND TO BE CONSIOBRED IN PNDB CALCULATION HORMALL? 1 

C 12 = LAST BAND 

C DB = ARRAY CONTAINING THE SPBCTHUM 9HOSE PNL IS TO BE CALCULAPED 

C PNDB1 = PNL CALCULAPED BY SMOOPHED CURVE FIT OF THE TABULAR 

C DATA 

C PNDB2 = PNL CALCULATED WSINS TABULAR DArA. 

C 

DIMENSION DB1 (9) ,DB10 (9) ,DB64 (15) ,DB50 (10) , DB(27) 

DATA (DB1 (I) ,I=1,9)/63. 6, 59.7 , 55 .2,53. , 50.5 , 48. ,45. 6,44. 5,43./ 

DATA (DBIO (I) ,1=1,9) /87.,85.,83. 1 , 80 . , 79. 2, 78. , 75. , 75. ,74./ 

DATA (DB64 (I) ,1=1,15) /1 12. ,111. , 109., 10 7. ,106., 105., 103. , 102., 101. 

1, 100., 100., 100., 100. ,100., 98./ 

DATA (DB50 (I) ,r=1,10)/91. ,89, , 3 7 . , 86. , 86. , 87., 88. , 9 1. , 94. , 97. / 

SUMN1=0.0 

SUMH2=0.0 

FMAX1=0.0 

FHAX2=0.0 

DO 50 I=T1,I2 

IF (I .GT. 15) SO TO 30 

EX=0 . 0301 03* ( DB ( I) - DB 64 ( I ) ) 

FN1=64.0*10.0**EX 

FN2=FN1 

IF (FN1 ,GT. 15. .OR. I . GT. 9) GO TO 40 
EX=(DB(I) -DB1 (I) ) /(DB10 (D -DB1 (I) ) 

FNLOH=10.0**EX 

IF (FNLOW .LT. FNl) FN1=FNLOW 
FN2=FN1 
GO TO 40 
30 J=I-15 

DBC=DB50(J) 

EX=0,030103* (DB (I) -D8C) 

FN1=50.0*10.0**EX 

IF (DB(I) .GE. DBC) GO TO 35 

X= (OBC-DB(I) ) **3 

FN1=FN1*(1.0+0.0279*X/(1.0+0.6*X) ) 

FN2=FN1 
GO TO 40 

3 5 X=DB(I)-DBC 

AEG=0.2*3. 1415927*X 
FS=SIN (ARG) 

COB=FS 

IF (X .GT. 10.0) GO TO 31 
C0R= (3.5-0.275*X) *FS 

IF (X .LT. 3.0) COR=COR+0. 6*X*(3.0-X) 

37 FN2=FN1* (1.0 + 0. 01*C08) 

40 SDHN1=SUMN1+FN1 
SDMN2=SUMN2+FN2 
IF (FNl .GT. FHAX1) FMAX1=FNl 
IF (FN2 .GT. FMAX2) FMAX2=FN2 
50 CONTINUE 

Fl=0.85*PMAX1+0. 15*SUMN1 
F2=0. 85*FMAX2+0. 15*SUMN2 
PNDB1=10.D*(4,0+ALOG10 (FI) /O. 30103) 

PNDB2 = 1 0.0* (4.0 + AL3G 10 (F2) /O. 30103) 

RETURN 

END 
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SIBFPC aET2 

SOBROOTINE JET (A , AH, V J, S JET) 

C 

r THIS SOBROOTINE CALCnLArBS JET NOISE 1/3 OCTAVE SPECTRA AT 50 AH5 120 

C DEGREES {SJET(I, 1) AND SJET(I,2) RESPECTIVELY) ON A 100 FT. SIDELISE. 

C CALLS BASPAT AND ASSOCTATED S03R00TINES. 

C ARGOMENTS 

C A = NOZZLE AREA, SQOAEE FEET. 

C AH = JET CHARACTERISTIC LENGTH, INCHES 

C VJ = JET EXIT VELOCITY, FT/SEC 

C SJET = JET SPECTRA, DB 

COHMON /FREQ/ F(27) 

COMMON /COND/ DKIND , AKIND, T, RH, D 1ST 

COMMON /CONSTS/ NF,RPM, ROTORB , FL OH, VR, A MACHF, AM ACHR, 3 OC 
COMMON /ANGLES/ ANG1,ANG2 
COMMON /SPECTR/ SJ(27,2) 

DIMENSION DI (16) ,ANG (16) 

DIMENSION SJET (27, 2) 

DIMENSION OAJET(2) 

DIMENSION B(5) 

DIMENSION ATA (27), R (2) 

DIMENSION TFA(27,2) 

DATA BK1 , BK2/. 05 , 1 . /, NDG/U/ 

DATA AO, A1/19.9,26.3/ 

DATA BO /-4.7515/, (B(I),T=1,4) /. 60049, -22. 942, 4. 3 15, 1 1. 1 46/ 

DATA CO,C1/-10.8,-15.82/ 

DATA (DI (I) ,1=1 ,16) /-7.4,-6. 8,-6. 3,-5. 7, -5. 2, -4. 4, -3. 6, -2. 7,-1 .8, 
1-1.,0.,.B, 1.8, 3.1, 4.5,6. /,REFDC /2.5/, REFANG /135./ 

DATA (ANG (I) ,1 = 1,15) /10. ,20. , 3 3 . , 40. , 50. , 6 0. , 70. , 80. , 90. , 1 00. , 1 1 0 
1.,120. ,130.,140.,150.,160./ 

PI=3. 1415927 
RT=200. 

NB=3 

ESTD=100. 

COMTRM=93.*ALOGlO(VJ) - 1 82 . 2+ 1 0. * AL031 0 (A) 

FIFTY=PI*ANG1/180. 

FORTY F=P I* R EFA N G/1 80. 

SIXTY=PI*ANG2/180. 

R (1)=DIST/SIN (FIFTY) 

R (2) -DIST/STN (SIXTY) 

DOWN1=ANG1-.01 

DOHN2=ANG2-.01 

0P1=ANG1+,01 

OP2=ANG2+.01 

DO 10 1=1,16 

I=I 

IF (ANG (I) .GT. DOWN1 .AND. AN3(I) . LT. DPI) KEEP1=I 
IF (ANG (I) .GT. DOWN2 .AND. ANG (I) . LT. 0P2) KEEP2=I 
10 CONTI NHE 

ADJST1=DI (KEEP1) -BEFDI 
ADJST2 = DI(KBEP2) -REFDI 

OAJET (1) =COMTRM + 20. *ALOG10 (2. *3IN (FIFTY) /SIN (FORTY?) ) ♦ADJSTI 
OAJET (2) =COMTRM + 20.*ALOG10 (2.*SIN (SIXTY) /SrN(FORTYF) ) 4-ADJST2 
THaDOC=10.*ALOGl0 (3.) 

PACTR= (AH/12.) /VJ 
DO 30 1=1, NF 
1=1 

STBHAL=F(I) *FACTR 
X=ALO310 (STEHAL) 
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IF (STBHSL .GT. BK1) GO TO 20 

DECB=AO+a1*X 

GO TO 25 

20 SaH=BO 

IF (STRHAL .37, BK2) GO TO 22 
X=X-ALOG10 (. 23) 

DO 21 IDG=1,NDG 
IDG=IDG 

21 SUM=SOM+B (IDG) * (X**IDG) 

DECR=SOa 

GO TO 25 

22 DECR=CO+C1*X 

25 SJET(I,1)-OAJET{1)+DECR-rHEDOC 
30 SJET (1, 2) =OAJET(2) +DRCa-THRDOG 
DO 40 J=1,2 

40 CALL BASPAT (SJET ( 1, J) ,RT ,NF, S3 , T, RH, TF A (1 , J) , B ( J) , AT A| 

COMBINED ADJTJSTHEMT (ATM ATTN) IN - TO 100 FOOT RADiaS AT 45 
DEGBEES, THEN OUT - TO 100 FOOT SIDELINE AT 50 AND 120 DEGREES. 

DO 50 1=1, NF 

SJET (I, 1) = SJET(I, 1) +TFA (1,1) T (RT/SIN (FDBTYF) -RSTD/SIN (FIFrY) ) /I 000 

1 . 

50 SJET (1,2) =SJET(I,2) +TFA (1,2) * (JT/SIN (FDRTYF) -RSTD/SIN (SIXTY) ) /1000 

1 . 

ADJUST JET SPECTRA TO DESIRED (DI3T) SIDELINE. 

DO 60 1=1, NF 

SJ (1,1) =3 JET (I, 1) +TEA (1,1) * ( (RSTD-DIST) /SIN (FIFTY) | /I 300. 

6 0 SJ(I,2) =S JET (1,2) tTFA (1,2) * ( ( RS T D-DIST) /SIS (SIXTY) ) /1 330. 

DELDIS=20.*ALOG10 (DIST/RSTO) 

DO 61 1=1, NF 
SJ (I , 1 ) = S J (1,1) -DELD IS 
6 1 SJ (1,2) =SJ (1,2) -DELDIS 

PRINT OUT JET NOISE SPECTRA. 

WRITE (6,101) 

101 FORMAT (1H0) 

WRITE (6,100) 

103 FORMAT (56X,17HJET NOISE SPECTRA//) 

WEITE (6,110) 

113 FORMAT (44X,43HAT 103 FT. SIDELINE, AND 50 AND 120 DEGREES//) 

CALL XPLANE 
WHITE (6,102) 

102 FORM AT (53X,30H FREQ SJETF SJETP //) 

DO 105 1=1, NF 

105 WRITE (6,106) F ( I) , S JET ( 1 , 1) , S JE T (I , 2) 

105 FORMAT (50X,F10. 0 ,2F10. 1) 

WRITE (6,1000) 

1030 FORMAT (1H1) 

RETURN 

END 
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SIBFtC BASPTT 

SOBHOOTINE BASPAT ( A, RT, NF, NB, P, RH, TFA, R, ATA) 

C CALLS ATMAT 

DIMENSION A (27) ,TFA (27) , ATA (27) 

REAL M(27,2) 

NLS=NF-1 

DO 10 1=1, NLS 

M(I,2)=A(I+1)-A(I) 

ID M (1+1, 1) =M (1,2) 

M(1,1)=M(1,2) 

M(NP,2)=M(NF,1) 

NS=3 
S = NS 
B=NB 

BI=B*2.0*FLOAT (NS) 

F 1 = 1 . 0 

IF(NB-EQ, 1) F1=10.0** 1.8 
IF (NB.EQ-3) F1=10.0**1.7 
DO 30 1=1, NF 

FC=10. 0** (O. 3/3*FLDAT (1-1) ) *F1 
CALL ATMAT (T, RH , 1 000 . O , FC, AC) 

SL=1-0 

S1A=10.0** (- AC/1 0.3) 

DO 20 K=1,NS 

XL=-K 

XR=K 

EL=M (I,1)/2.0*XL/S 
ER=M (I, 2) /2.0*XR/S 

SL=SL+10.0** (EL/10.0) 4^1 0. 0** (ER/10. 0) 

FL= 1 0. 0** ( 0 . 3/B I* XL ) * FC 
FR=10.0**(0.3/B1*XR) *FC 
CALL ATMAT (T ,R H, 1 000. 0 , FL, AL) 

CALL ATMAT(T,RH, 1 000. 0, FR , AR) 

20 S1A=SLA+10.0** ( (EL-AL)/10.0) +13. 0** ( (ER-ARl /10.0) 
TFA(I) =10.0*AL3G10(SL) - 1 0. 0* ALD3 1 0 (SLA) 

30 ATA{I)=A (I)- TFA(I) *R/1 000. 0-20. 0*ALOG1 0 (R/RT) 
RETDRN 
END 


$lB?rC ATMATV 

SOBROOTINE ATMAT(T,RH,DIST,FRE3, ATT) 

DIMENSION A{22) 

DATA A/0.870, 0.750,0.652,0.570,0.505,0. 452,0.406,0.369,0.335, 

1 0.308,0.286,0.268,0.253,0.240,0.231,0.225,0.220,0.215,0.210, 

1 0.208,0.202,0.200/ 

AC=(0. 1* (PR B3/1 000.0) **2.05) /( 1.651 -.001 03*T) **2.05 
AMM= (10.0* (FRBQ/1 000.0) **1.003) /10. 0** (0. 52-. 00 504* (T + 

1 SQRT(256.0- (10.0-T/5.0)**2) ) ) 

HA=0.25*RH/10.0**(1. 493-.01638*T-.02*SQRT (128.2- (1 0. O-T/5. 00 ) **2) ) 
HMM=10.0** (0.4973*ALOG10 (FREQ) - 1 . 4894) 

HH=HA/HMM 

IF (8H.GT.0.25) GO TO 10 
AA=1.2*HH 
GO TO 15 

ID IF(HH.GT.0.60) GO TO 
AA=1.543*HH-.086 
GO TO 15 
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11 IF (HH.GT.0.95)SO TO 12 

AA=0.84+0. 16*SIS ( 3. 1 4159/2. 0* (39-0.6) /0.35| 
GO TO 15 

12 IF(HH.GT. 1.25) GO TO 13 

AA=0.87 + 0. 1 3*COS (3. 141 59/2. 0* (HH-0. 95) /O. 3) 
GO TO 15 

13 IF (HH.GT.6.5) GO TO 14 
HTE3T=1.25 

DO 16 1=2,22 

HTEST=HTEST+0.25 

IF (HH.LE.HTEST) GO TO 17 

16 CONTINUE 

1 7 AA=A (I) + ( (HTEST-HH) /O. 25) * (A (1-1 ) -A (I) ) 

GO TO 15 

14 A A =0.2 

15 CONTINUE 

ATT= (AMfl*AA+AC) * (DIST/1 000 . 0 ) 

RETURN 

END 


SIBFTC SIDLTT 

SUBROUTINE SIDLAT (NM, AI , RSTD, S I DIST, SIN AT, B ADIST) 
DIMENSION A1 (19) ,SINAT (19) ,RADtST (19) 

F=3. 1415927/180.0 

eSQ=2 0.0’» ALOG10 (SIDIST/RSTD) 

DO 10 J=1,NH 
ST=SIN (AI (J) *F) 

IF(ST.LE.O.O) GO TO 20 
SINAT ( J) =RSQ-20. 0*ALOS10 (ST) 

RADIST (J) =SIDIST/ST 
SO TO 10 
20 SINAT (J) =0.0 
SADIST (J) =0.0 
10 CONTINUE 
RETURN 
END 


SIBFTC XPLANY 

SUBROUTINE XPLANB 

SIMPLE PRINT-OUT ROUTINE 

COMMON /XPLAN/ XPLN1 (15) , XPLN2 ( 1 5) 
WRITE (6,3) (XPLNI(I) ,1=1,12) 

3 FORMAT (45X,12A6//) 

WRITE (6,3) (XPLN2(I) ,1=1,12) 

RETURN 

END 
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SUBROUTINE &TSG 
PURPOSE 

NDIN POINTS A GIVEN GENERAL TABLE ARB SSLECrED AND 
ORDERED SUCH THAr ABS {ARG ( I) -X) . 3E« ABS (ARG { J|--X) IF I. ST. !• 


USAGE 

CALL ATSG (X, Z /WORK, I ROtf , ICOL, ARG, V AL, S DIE) 


DESCRIPTION 

X 

z 

p 


WORK 
I ROW 

I COL 
ARG 

VAL 


NDIM 


OF PARAHETER5 

THE SEARCH AR3UEENT. 

THE VECTOR OF ARGUMENT VALUES {DIHENSIDN XROW). 

IN CASE ICOL=1, P IS THE VECTOR OF FUNCTION VALUES 
(DIMENSION IROW) . 

IN CASE ICOL=2, F IS AN IROW BY 2 MATRIX. THE FIRST 
COLUMN SPECIFIES THE VECTOR OF FUNCTION VALUES AND 
THE SECOND THE VECTOR OF DERIVATIVES. 

A WORKING STORAGE (DIMENSION I ROW). 

THE DIMENSION OF VECTORS 2 AND WORK AND OF EACH 
COLUMN IN MATRIX F. 

THE NUMBER OF COLUMNS IN F (I.E. 1 OR 2). 

THE RESULTING VECTOR OF SELECTED AND ORDEREO 
ARGUMENT VALUES (DIMENSION NDIM) • 

THS RESULTING VSCTOR OF SELECTED PUSCriOH VALUES 
(DIMENSION NDIM) IN CASE IC0L=1. IN CASE ICOL=2, 

VAL IS THE VECTOR OF FUNCTION AND DERIVAriVE VALUES 
(DIMENSION 2*NDIM) WHICH ARE STORED IS PAIRS (I.S. 
EACH FUNGTIDN VALUE IS FOLLOWED BY ITS DERIVATIVE 
VALUE) • 

THE NUMBER OF POINTS WHICH MOST BE SELECTED OUT OF 
THE GIVEN TABLE (2,F). 


REMARKS 

NO ACTION IN CASS IROW LESS THAN 1. 

IF INPUT VALUE NDIM IS GREATER THAN IROW, THE PROGRAM 
SELECTS ONLY A MAXIMUM TABLE OF IROW POINTS. rflEBEFORE THE 
USER OUGHT TO CHECK CORRESPONDENCE BETWEEN TABLE (ARG, VAL) 
AND ITS DIMENSION BY COMPARISON OF NDIM AND IROW, IN ORDER 
TO GET CORRECT RESULTS IN FORTHEa WORK WITH TABLE (ARG, VAL). 
THIS TEST MAY BE DONE BEFORE OR AFTER CALLING 
SUBROUTINE ATSG. 

SUBROUTINE ATSG SSPECIALL? CAN BE USED FOR GENERATING THE 
TABLE (ARG, VAL) NEEDED IN SUBROUTINES ALI, AHI, AND ACPI. 


C SUBROOTINES AND FUNCTION SUBPROGRAMS REQUIRED 

C NONE 

C 

C METHOD 

C SELECTION IS DONE BY GENERATING THE VECTOR WORK WITH 

C COMPONENTS WORK (I) = AB3 ( Z ( I) «X) AND AT EACH OF THE NDIM STEPS 

C (OR IROW STEPS IF NDIM IS GREATER THAN IROW) 

C SEARCHING FOR THE SUBSCRIPT OF THE SMALLEST COMPONENT, WHICH 

C IS AFTERWARDS REPLACED BY A NUMBER GREATER THAN 

C MAX (WORK (I) ) . 
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SOBROOriNE &TSG(X, 3,F,WORK,iaOI, I30L,lRS,V!kL,NDtH( 

THIS SUBROUTINE HAS BEEN ALTERED TO HANDLE PATHOLOGICAL SITUATIONS SUCH 
AS EXTRAPOLATIONS WITHIN, OR OUTSIDE, OF THE TABLE. IF EXTRAPOLATION 
OUTSIDE THE TABLE IS REQUIRED, NOIM ON RETURN IS REDUCED TO 2 AND A 
STRAIGHT LINE THROUGH THE TWO CLOSEST POINTS IS USED TO DO THE 
EXTRAPOLATION. OTHERWISE ARG IS FORGED TO BRACKET THE VALUE OF X. 


DIMENSION Z(1),P(1) , WORK (1), 183(1), VAL(1) 

NSAVE=NDIM 

NDIM=IROW 

IP (IRON) 11 ,11,1 

1 N=NDIM 

IF N IS GREATER THAN IRON, N 13 SET EQUAL TO IRDN. 

IF (N-IROH) 3,3,2 

2 N*IROH 

GENERATION OF VECTOR WORK AND COMPOTATIOH OF IIS GREArESr ELENENr. 

3 B=0. 

DO 5 1=1, IRON 
DELTA=ABS (Z (I) -X) 

IF (DELTA-B) 5,5,4 

4 B=DELrA 

5 WORK {I) =DELTA 

GENERATION OF TABLE {ARG,VAL) 

B=B+1. 

DO 10 J»1,N 

DELTA=B 

DO 7 1=1, I ROW 

IF (WORK (I) -DELTA) 6,7,7 

6 II=I 

DELTA=WOHK (I) 

7 CO NT IN HE 
ARG (J) =2 (II) 

IF(ICOL-I) 8,9,8 

B VAL(2*J-1)=F(TI) 

III=II+IROH 
VAL {2*J) =F (III) 

GOTO 10 
9 VAL(J)=P(I1) 

10 WORK (II) =B 

DO CALCULATION TO DETERMINE WHETHER X, THE SEARCH ARGUMENT, IS OUTSIDE 
THE RANGE OF THE TABLE. IF SO, EXTRAPOLATION IS REQUIRED AND THE 
SIMPLEST TYPE OF CURVE, A STRAIGHT LINE, SHOULD BE USED IN ORDER TO 
AVOID THE LIKELIHOOD OF UNREASONABLE RESULTS FROM A POLYNOMIAL OF 
HIGHER DEGREE. 

ILT=1 

IGT=1 

DO 15 I=1,IR0W 
IF (X .GT. Z(D) ILT=0 
IF (X .LT. Z(D) IGT=0 
15 CONTINUE 
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NOW, IF X IS GREATER THAN ANY ONE OF THE TABLE ARGUMENTS, ILT=0 
NOW, IF X IS LESS THAN ANY ONE OF THE TABLE ARGUMENTS, TGT=0 
HOWEVER, IF ILT REMAINS T, THEN X IS LESS THAN ALL ARGUMENTS OF THE TABLE 
IF IGT REMAINS 1, THEN X IS GREATER THAN ALL THE ARGUMENTS OF THE TABLE. 
IN EITHER OF THE LATTER EVENTS, THE TWO CLOSEST VALUES FROM THE TABLE 
SHOULD BE USED FOR THE EXTRAPOLATION, SINCE ARG IS ORDERED SUCH AS TO 
BE AT INCREASING DISTANCES FROM X, ONE MERELY HAS TO LIMIT CONSIDERATION 
TO THE FIRST TWO VALUES OF ARG AND VAL. 

IF (ILT. EQ. 1) GO TO 17 
IF (IGT .EQ. 1) GO TO 17 
GO TO 18 

17 NDIM=2 
RETURN 

CHECK TO MAKE SURE THAT THE RANGE OF ARG BRACKETS X FOR THE CASE WHERE 
X FALLS WITHIN THE RANGE OF THE TABLE OF ARGUMENTS, 2(1). 

THIS IS TO AVOID UNNECESSARY EXTRAPOLATION WITHIN THE TABLE FOR CERTAIN 
SPECIAL CASES OF THE SPACING OF Z(I). 

18 CONTINUE 
ILT=1 
IGT=1 

DO 21 I=1,NSAVE 
IF (X .GT. ARG(I)) ILT=0 
IF (X .LT. ARG(I)) IGT=0 
21 CONTINUE 

IF (ILT .EQ. 1) GO TO 23 
IF (IGT .EQ. 1) GO TO @( 

NDIM=NSAVE 

RETURN 

23 NSP1=NSAVE+1 
DO 25 I=NSPl,IROW 
IF (X .GT. ARG(I)) GO TO 26 

25 CONTINUE 
RETURN 

26 IWHICH=I 
ARG(NSAVE)=AR6(IWHICH) 

VAL(NSAVE)=VAL(IWHICH) 

NDIM=NSAVE 

RETURN 

29 NSP1=NSAVE+1 

DO 33 I=NSP1 ,IROW 

IF (X .LT. ARG(I)) GO TO 34 

33 CONTINUE 
RETURN 

34 IWHICH=I 
ARG(NSAVE)=ARG(IWHICH) 

VAL(NSAVE)=VAL(IWHICH) 

NDIM=NSAVE 

11 RErURN 
END 
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SUBROUTINE LGRTRP (XX,YY.UM,X,Y) 


INTERPOLATE A FUNCTION VALUE Y AT DESIRED ARKJMENT VALUE X FROM A TABLE 
OF VALUES YY (TABULATED Y) VERSUS XX (TABULATED X). 

XX(I) MUST BE ARRANGED IN ORDER OF INCREASING DISTANCE FROM X AS I 
INCREASES. 

XX MUST BRACKET X, UNLESS EXTRAPOLATION IS CALLED FOR. 

MODIFIED ATSG CAN BE USED TO ACCOMPLISH THESE ORGANIZATIONS PRIOR TO 
calling LGRTRP. 


DEFINITION OF PARAMETERS IN SUBROUTINE CALL 
XX ARRAY CONTAINING TABULATED ARGUMENT (X) VALUES 

YY ARRAY CONTAINING TABULATED FUNCTION (Y) VALUES 

LIM NUMBER OF XX VALUES IN THE TABLE 
X VALUE OF ARGUMENT AT WHICH FUNCTION VALUE IS DESIRED 
Y CALCULATED DESIRED FUNCTION VALUE. 

METHOD IS LAGRANGE INTERPOLATION AS GIVEN ON P 2T5 OF NUMERICAL 
METHODS AND COMPUTERS BY S.S. KUO, PUBLISHED BY ADDISON-WESLEY, 1965 

DIMENSION YY(12),XX(12) 

NDIM=LIM 
3 CONTINUE 
C=0. 

DO 52 1=1 ,LIM 
ZXI=XX(I) 

P=l. 

DO 40 J=1,LIM 
IF (I-J) 35,40,35 
35 ZXJ=XX(J) 

A=(X-ZXJ)/(ZxI-ZXJ) 

P=P*A 

40 CONTINUE 
B=P*YY(I) 

C=C*B 

52 CONTINUE 
Y=C 


CHECK FOR RIDICULOUS RESULTS IN INTERPOLATION, USE STRAIGHT LINE 
INTERPOLATION IF HIGHER ORDER INTERPOLATION GIVES VALUE OF Y WHICH IS 
GREATER OR LESS THAN BOTH THE YYS AT THE NEAREST TWO XXS. NOTE IF THE 
FUNCTION IS KNOWN NOT TO BE MONATONIC, THEN THIS LOGIC SHOULD PROBABLY 
BE DISABLED. 

IF (NDIM .EQ. 2) GO TO 82 
DELTA=X-XX(1 ) 

DO 72 1=2, LIM 
DELTA2=X-XX(I) 

IF (DELTAl .LT. 0. .AND. DELTA2 .GT. 0.) GO TO 73 
IF (DELTAl .GT. 0. .AND. DELTA2 .LT. 0.) GO TO 73 

72 CONTINUE 

73 IWHICH”! 

IF (Y .GT. YY(1) .AND. Y .GT. YY(IWHICH)) GO TO 75 
IF (Y .LT. YY(1) .AND. Y .LT. YY(IWHICH)) GO TO 75 
GO TO 82 
75 LTM=2 

XX (2)=XX(IWHICH) 

YY (2)=YY(IWHICH) 

GO TO 3 
82 CONTINUE 
LIM=NDIM 
RETURN 
END 
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Printed Output 


The variable names used in the program output have been selected for clarity and 
ease of interpretation. In general the trailing letters ’*F" and ”R” of a given name 
(e. g. , STONEF) refer to the sound level at the forward and rear angles of maximum 
noise, respectively. Furthermore, the leading *'S” stands for ’’sound”. Thus 
STONEF is the name of the sound spectrum of tones at the forward and maximum angle. 
Otherwise symbolism should be self-explanatory. For example PR stands for pres- 
sure ratio and ATTN stands for attenuation. 

The initial page of the output section is a set of input data in proper format to be 
used in running the program. 


TABLE Vn. - SET OF INPUT DATA-NOISE ESTIMATE 

INPUT oata 

HAHLTCN STAMDAI^C 1.2 PP PPOP-FAN (QF-9) 

100 PEPCENT SP5EC, TAKEOFF BL AOF-S E TT I NG 


890.0 700.0 2.0 

579.00 21. 76 18.60 

S 

500.0 59.0 70.0 

25.0 25.0 

90.00 


2227.0 


15.0 


12.5 


TABLE Vm. - CALCULATED BLADE PASSING TONE AND ITS HARMONICS 

AT 30. 5~METER (100-FT) SIDELINE AND 50° AND 120° 

[Hamilton standard 1.2 PR prop-fan (QF-9); 100 percent speed; takeoff blade 
setting. ] 


f P^Q 

S^roNEF 

STCl^rp 

5 00. 

96.4 

0 2.8 

1000. 

90.4 

96. B 

16 00. 

86.9 

93.3 

20 00. 

64. 4 

90.8 

2500. 

82. 5 

8 8.9 

3150. 

60.9 

8 7.3 

40 00. 

82.0 

88.4 

5000. 

79.9 

86.3 

6 300 . 

78.3 

84.7 

80 00. 

79.3 

85.7 

10000. 

77.2 

83.6 
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TABLE DC. - FANBROADBAND SPECTRUM AT 30.5-METER (100- FT) 
SIDELINE AND 50° AND 120° 

[Hamilton standard 1.2 PR prop-fan (QF-9); 100 percent speed; takeoff 
blade setting . ] 


FR5Q 

S0ROAO 


5^r 

63* 

59.0 

80. 

62.9 

ICO. 

66 . 8 

125. 

71.0 

160. 

75.6 

200. 

79 . 7 

250. 

83.6 

315. 

87.1 

400. 

90.0 

500. 

91.8 

630. 

92.8 

.800. 

-92it8, 

1000. 

91.8 

1250. 

89.9 

1600. 

86.9 

2000. 

83.6 

2 500. 

79.7 

3150. 

75.5 

4000. 

71.0 

5000. 

66.8 

6 300. 

62.8 

80 30. 

58.9 

lOOOO. 

55.7 


TABLE X. - JET NOISE AT 30. 5-METER (100- FT) SIDELINE 
AND 50° AND 120° 

[Hamilton standard 1. 2 PR prop-fan (Q F-9); 100 percent speed; 
takeoff blade setting. ] 



SJETF 

SJETR 

50. 

76.2 

83.2 

63. 

77.1 

84.2 

80. 

77.6 

84.6 

100. 

77.6 

84.6 

125. 

77.1 

84.2 

160. 

76.3 

83. 3 

200. 

75.2 

8 2.3 

253. 

73.9 

81.0 

315. 

72.6 

79.7 

400. 

71.2 

78.2 

500. 

69.7 

76.7 

6 30. 

68.1 

75.2 

800. 

66.5 

73.6 

1000. 

65.0 

7 2.1 

1250. 

63.5 

70.6 

1600. 

61.9 

69.0 

2000. 

60.5 

67.6 

2500. 

59.1 

66.2 

3150. 

57.7 

64.9 

4000. 

56.4 

63.6 

5000. 

55.4 

62.6 

6300. 

54.5 

61.8 

8000. 

53.9 

61.3 

10003. 

53.8 

61.4 



TABLE XI. - TOTAL SPECTRA AT 30.5-METER (100- FT) SIDELINE 


AND 50° AND 120° 


[Hamilton standard 1.2 PR prop-fan (QF-9); 100 percent speed; 
takeoff blade setting; internally generated forward OASPL, 
101.5; internally generated rear OASPL, 105.7.] 


FREQ 

5PLF 

SPLF 

50. 

76.2 

83.3 

63. 

77.2 

84. 2 

80. 

77.7 

84.7 

100. 

77.9 

84.7 

1 25. 

78. 1 

84.4 

160. 

79.0 

84.0 

200. 

81.1 

84.2 

2 50. 

84.0 

85.5 

315. 

87.2 

87.8 

400 . 

90.0 

90.2 

500. 

96.5 

102.9 

6 30. 

92.9 

9 2.9 

800. 

92*8 

92.9 

1000. 

91.8 

96.8 

12 50. 

90 . 0 

9C. 0 

1600. 

87.0 

93*3 

2000. 

84.4 

90. 8 

2500. 

82.5 

88.9 

3150. 

80.9 

87.3 

4000. 

B2.0 

88*4 

5000. 

79.9 

86*4 

6300. 

78.3 

84.7 

8000. 

79*3 

3 5.7 

100 00. 

77.2 

83.6 


TABLE Xn. - ATMOSPHERIC ATTENUATION PER 305 METERS 
(1000 FT) FOR SPECTRA OF TABLE X (ACCOUNTING FOR 
SPECTRUM SHAPE, FRONT AND REAR) 


FREQ 

DEL DBF 

D6LDBR 

50. 

0.1 

0.1 

63. 

0.1 

0.1 

80. 

0.1 

0.1 

100. 

0.1 

0.1 

125. 

0.2 

0.2 

160. 

0.2 

0.2 

200. 

0.3 

0.3 

250. 

0.4 

0.4 

315. 

0.5 

0.5 

400. 

0.6 

0.6 

500. 

0.7 

0.7 

630. 

0.9 

0.9 

800. 

1.2 

1.2 

1000. 

1.5 

1.5 

1250. 

1.9 

1.9 

1600. 

2.3 

2.4 

2000. 

3.1 

3.1 

2500. 

4.1 

4*1 

3150. 

5.6 

5.6 

4000. 

7.7 

7.7 

5000. 

io.8 

10.8 

6300. 

15.6 

15.6 

8000. 

22.4 

22.4 

100 00. 

31. 7 

31.7 



TABLE Xra. - TOTAL SPECTRA FOR TARGET PNL AT 

152. 5- METER (500- FT) SIDELINE AND 50° AND 120° 

[Hamilton standard 1.2 PR prop- fan (QF-9); 100 percent 
speed; takeoff blade setting; distance condition = side; 
atmospheric condition = atm. ] 


f p£g 

SPtF 

SPLP 

50* 

62.2 

69.2 

63. 

63.1 

70.2 

80. 

63.7 

7 j. 6 

100. 

63.9 

70.7 

125. 

64.0 

70.3 

160# 

64. 9 

69.9 

200. 

66.9 

70.1 

250. 

69.9 

71.3 

315. 

73.0 

73.6 

400. 

75.7 

76.0 

500. 

8 2.1 

88.5 

630. 

78.4 

78.5 

8 00. 

78.3 

78.4 

1000. 

77. 1 

82.2 

1250. 

75. G 

75.2 

1600. 

71.8 

78.2 

2000. 

68.8 

75.4 

2500. 

66.4 

7 3.0 

3150. 

64.0 

70.7 

40 00. 

64.0 

70.9 

50 00. 

60. 3 

67.4 

6300. 

56.2 

63.5 

80 00. 

53.6 

61.4 

10000. 

46.7 

55.0 


TABLE XIV. ~ HARD WALL PNL'S AT NOTED POINT 

[Hamilton standard 1.2 PR prop-fan (QF-9); 100 per- 
cent speed; takeoff blade setting. ] 


PN08F = 

95.2 

PNOBR = 

mo. 5 

300. 

SIDE 
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table XV. - ATMOSPHERIC ATTENUATION SPECTRA PER 
305 METERS (1000 FT) FOR JET NOISE SPECTRUM 
(ACCOUNTING FOR SPECTRUM SHAPE, FRONT AND REAR) 


fRHQ 

50. 
63. 
80. 
100 . 
125. 
160. 
200 . 
250. 
315. 
400. 
500. 
6 30. 
800. 
1000. 
12 50. 
1600. 
2000 . 
2500. 
3150. 
40 30. 
50 00. 
6300. 
8000. 
13000. 


DEL DBF 

0.1 
0 . 1 
0.1 
0.1 
0.2 
0*2 
0.3 
0.4 
0.5 
0.6 
0.7 
0-9 
1.2 

1.5 

1.9 

2.4 
3.1 
4 . 1 

5.6 

7.7 

10.9 

15.5 

22.4 
32,2 


DELOBP 

0.1 

O.l 

0.1 

0.1 

0.2 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.9 

1.2 

1.5 

1.9 

2.4 
3.1 

4^1 

5.6 
/.7 

10.9 

15.5 

22.4 

32.3 


TABLE XVI. - JET SPECTRA FOR TARGET PNL AT 152. 5- 
METER (500- FT) SIDELINE AND 50° AND 120° 


[Hamilton standard 1. 2 PR prop-fan (QF-9); 100 percent 
speed; takeoff blade setting; distanee condition = side, 
atmospheric condition = atm. ] 


FREQ 

SPLF 

SPLP 

50. 

62.2 

69.2 

63. 

63.1 

70.2 

80. 

63.5 

70.6 

100. 

63.5 

70.6 

125. 

63.1 

70.1 

160. 

62.2 

69.3 

200. 

61.1 

68. 1 

250. 

59.8 

66.9 

315. 

58.4 

65.5 

400. 

56.9 

64. 0 

500. 

55.3 

62.4 

630. 

53.6 

60.8 

800. 

51.9 

59. 1 

1-000. 

50.3 

57.4 

1250. 

48.6 

55.8 

1600. 

46.7 

53.9 

2000. 

44.9 

52.2 

2500. 

43.0 

50.4 

3150. 

40.9 

48.3 

4000. 

38. 4 

46.1 

5000. 

35. 7 

43*6 

6300^ 

32.4 

•tO.6 

8000. 

28.2 

36.9 

10000. 

23.0 

32.5 



TABLE XVn, - CONSTANT NOY SPECTRUM REQUIRED 
TO REACH PNL GOAL (90.0 PNdB) 


[Constant Noy = 7.2.] 


FREO 

SGOAL 

50. 

81.2 

6 3. 

79.6 

80. 

77.9 

iOO. 

16.3 

12 5. 

74.6 

160. 

13.1 

200. 

71.7 

250. 

70. 3 

315. 

68.9 

400. 

68.5 

500. 

68.5 

630. 

68 . 5 

800. 

68.4 

1000. 

67.7 

1250. 

65.9 

1600. 

63.2 

2000. 

6Q»4 

2500. 

58.5 

3150. 

57.6 

4000. 

57.6 

5000. 

58.5 

6300. 

60.4 

8000. 

63.2 

10000. 

65.9 



TABLE XVm, - ATTENUATION SPECTRA REQUIRED TO 
REACH GOAL (90.0 PNdB) TOGETHER WITH ESTIMATED 
IN-DUCT OASPL'S AT 152. 5-METER (500- FT) SIDELINE 


[Hamilton standard 1. 2 PR prop-fan (QF-9); 100 percent 
speed; takeoff blade setting; in- duct forward OASPL, 
140.8; in-duct rear OASPL, 144.0.] 


FREQ 

SATTNF 

SATTNF 

50. 

-19.0 

-12.0 

63. 

-16.4 

-9.4 

80. 

-14.2 

-7.3 

100. 

-12.4 

-5.6 

125. 

“10.6 

-4.3 

160. 

-8.2 

— 3 . 1 

200. 

-4. 8 

-1.6 

2 50. 

-0.4 

1.1 

315. 

4. 1 

4.7 

400. 

7.2 

7.5 

500. 

13.6 

20.1 

630. 

9.9 

10.1 

800. 

9.9 

10.0 

1000. 

9.4 

14.5 

12 50. 

9.1 

9.2 

1600. 

8. 6 

15.1 

2000. 

8.4 

15.0 

2500. 

7.9 

14.5 

3150. 

6.4 

13.2 

4000. 

6.4 

13.3 

50 00. 

1.8 

8.9 

6300. 

-4.3 

3.1 

8000. 

-9.5 

-1.8 

10000. 

-19.3 

-11.0 

IN^OUCT 

nASPLF = 
OASPLP ^ 

140. 

144. 



APPENDIX B 


LINER SPECIFICATION PROGRAM DESCRIPTION 
Description 

This appendix presents a listing of the FORTRAN IV main program and subroutines 
required to specify the fan acoustic liners, once the required attenuation spectra have 
been determined. The variables calculated are tuning frequencies, lengths of liners, 
backing depths, and open area ratios. In the process decisions must be made for de- 
termining numbers of splitters, passage heights, facing sheet hole diameter, and facing 
sheet thickness. These latter variables are constrained by practical considerations, 
but some freedom is possible in iterating with the former list of variables to get a 
'’best'* result. This iteration requires some judgment and experience in order to be 
effective. 

The specific input data requirements are itemized at the beginning of the main pro- 
gram, and they will not be discussed in detail here. 

The printed output resulting from this program is presented in tables XK to XXni 
which appear after the main program and subroutines at the end of this appendix. 


Main Program and Subroutines 


c 

C SPECIFY LINER TUNING. L/H. AND ETA. 

C ASSUME REQUIRED ATTENUATION SPECTRUM IS GIVEN. 

C ACCEPTS A GENERAL NUMBER OF TUNINGS. UP TO 5. 

C CALCULATE OPEN AREA PATIOS AND BACKING DEPTHS* 

C 

C SUBROUTINES REQUIRED 

C ATTN 

C OFFOPT 

C BNOW 

C ANEWT 

C CR I T 12 

C ZRTDSl 

C ARBD 

C PNOBER 

C ATSG - MODIFIED IBM SYSTEM SUBROUTINE 

C LGRTRP - LAGRANGIAN INTERPOLATION 

C 
C 

C DATA REQUIRED 

C XPLN(I) 12A6 

C CONFIGURATION EXPLANATION CARD 

C PART, ETAOO 2F10.0 

C PART = MULTIPLE OF OPTIMUM THETA DESIRED FOR OFF-OPTIMUM DESIGN 

C ETAOO = ETA BELOW WHICH USE OFF-OPT IMUM DESIGN CONS IDERATIONS. 
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NFSWCHf NLINER 213 

NFSWCH... IF ZERO, TUNE FREQUENCIES ARE EXTERNALLY SPECIFIED 
IF NON^ZEROt FREQUENCIES ARE CALCULATED BY PROGRAM^ 

NLINER == NUMBER OF LINER TUNINGS DESIRED 
AM, T 2F10.0 

AM = DUCT MACH NUMBER 
T = TEMPERATURE IN DUCT. DEGREES R# 

FP(I) 5F10.0 

EACH VALUE OF TUNING GUESSED FOR ITERATION 

BLOHCD 5F10.0 

EACH VALUE OF L OVER H GUESSED FOR EACH TUNING 
H( I) 5F10#0 

DUCT HEIGHTS FOR EACH TUNING. IN. 

NOTE THAT YOU CAN USE PUNCHED OUTPUT FROM THE NOISE ESTIMATE 

PROGRAM FOR DATA BETWEEN ASTERISKS. OR GENERATE THESE DATA INDEPENDENTLY 

ISTART, ISTOP 213 

I START = FIRST NON-ZERO 1/3 OCTAVE BAND IN ATTENUAT ION GOAL SPECTRUM 
ISTOP = LAST NON- ZERO 1/3 OCTAVE BAND IN ATTENUATION GOAL SPECTRUM 
GOAL! I > 8F10.0 

GOAL = 1/3 OCTAVE SPECTRUM OF GOAL ATTENUATION. 

NF 13 

NUMBER OF 1/3 OCTAVE BANDS IN SPECTRUM 
IF ZERO. USED AS SWITCH. DO NOT READ SPL. SLJET. COND 
SPL(I) 8F10.0 

1/3 OCTAVE SPECTRUM OF NOISE GENERATION, FOR PNL CALCULATION 
SLJET(I) 8FI0.0 

1/3 OCTAVE spectrum OF JET NOISE, FLOOR FOR NOISE ATTN, FOR PNL CALC 
CONDin 12A6 

EXPLAINS LOCATION OF ABOVE SPECTRA 
OASPL FIO.O 

OASPL = DUCT OASPL. DB 

IF OASPL IS ZERO, IT ACTS AS A SWITCH, NO FURTHER READING AND NO 
BACKING DEPTHS OR OPEN AREAS ARE CALCULATED. 


PS FIO.O 

PS = DUCT STATIC PRESSURE, PSI 
THAin 5F10.0 

THA(l) - FACE PLATE THICKNESS FOR FIRST TUNING, IN. 
ETC 

DA(I) 5F10.0 

DAIl) = HOLE DIAMETER FOR FIRST TUNING, IN. 

ETC. 

COMMON /FREQ/ F(27) 

COMMON /PAR/ PART,ETAOO 

COMMON /PESANT/ PS.AM.T.OASPL.TH.D 

COMMON /FLOBY/ NBYSW 

DIMENSION DOB(27) ,GOAL(27) ,FU( 27» 

DIMENSION XPLNIIS) 

DIMENSION DIFFEIS ),ET(5) ,OIFFL (5) , BL0HI5) ,DDBN(5 ,27) 

01 MENS ION FP( 5) .H(5) .ETG( 5) .FPG< 5) .BLOHGIS ) . BL (5 ) . HG( 5 ) 
DIMENSION THETA(5),CHI(5),SIGMA(5),B(5) 

DIMENSION THA(5),0A(5) 
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DIMENSION SPU27) .SLJET<27).C0N0a5) 

DIMENSION PEAKIS) iPERLIS) 

NAMELIST /NMl/ NF ,FU, AM,T,C.ET»FP ,BLOH, I START, ISTOP, GOAL ,H,NFSWCH, 
INL INER.NLIN.PART. ETAOO 

NAMELIST ANM2/ AM ,ET ,FP,NF ,ODBN,DOB,FRAC 
NAMELIST /NM3/ IT ERAT, FR AC , TE STE , DI FFE ,TESTL ,D IFFL 
NAMELIST /NM4/ AM,T.PS»QASPL,THA,OA 
NAMELIST /NM5/ NF,SPL,SLJET 

NAMELIST /NM6/ AM ,C, PART , ETAOO .FRC .FP, ET .BLOW. PERL .PEAK. 0D8.D0BPK. 
IDDBN 

COMPLEX FIMPfEIG 

DATA ( FUdl .1 = 1.2 7) /50. .63. . 80. . 100. . 125. .1 60 .. 200. . 250. . 315. .400 

1.. 500., 630., 800., 1000., 1250., 1600. ,2000. ,2500. ,3150. ,4000. ,5000. ,6 

2300. . 8000. . 10000. .12500. . 16000. .20000. / 

TP, = 1. 

NBYSW=0 

NBYSW IS A SWITCH FOR CONTROL OF FLOW BY EFFECT IN IMPEDANCE MODEL 
CALCULATION 

= 0 USE ADDITIVE 

= 1 USE multiplicative 

6 CONTINUE 
NF=24 
NFO=NF 

DO 1 1 = 1.27 

1 F(I)=FU(I) 

1090 CONTINUE 

READ (5,100) (XPLNd ) .1=1.12) 

100 FORMAT (12A6) 

WRITE (6,999) 

WRITE (6.7) 

7 FORMAT (60X,10HINPUT DATA///) 

WRITE 16,1005) (XPLNII ),I=1,12) 

1005 FORMAT (35X.12A6//) 

READ (5,102) PART, ETAOO 
102 FORMAT (2F10.0) 

WRITE (6.1020) PART, ETAOO 
1020 FORMAT ( 35X , 2F10. 2// ) 

FRC=1./PART 

READ (5,101) NFSWCH, NLINER 

101 FORMAT (213) 

WRITE (6.1010) NFSWCH, NLINER 
1010 FORMAT {35X,2I3//) 

NLIN=NLINER 
READ (5.2) AM.T 

2 FORMAT (2F10.0) 

WRITE (6.200) AM.T 

200 FORMAT (35X ,F10.2 .FIO. 1// ) 

C=49.a2*SQRT(T) 

READ (5.3) (FP(I) ,I»1.NLIN) 

3 FORMAT (8F10.0)- 

WRITE (6,300) (FP(I),I*1,NLIN) 

300 FORMAT (35X.5F10.0) 

READ (5,3) (BLOH( I), I=1,NLIN) 

WRITE (6,301) (BL0H(I),I=1,NLIN) 

301 FORMAT (1H0.34X.5F10.1) 

READ (5,3) (H(I ), I=1,NLIN) 

WRITE (6,302) (H( I), I=1,NLIN) 

302 FORMAT ( 1H0,34X.5F10.2) 

READ (5,4) ISTART, ISTOP 
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4 FORMAT C2I3) 

WRITE (6.1011) ISTART,ISTOP 

READ (5,5) (GOAU I),I=ISTART,ISTOP) 

5 FORMAT (SFIO.O) 

WRITE (6,500) (GOAL( I ) ,I=ISTART,ISTOP) 
50D FORMAT ( IHO , 34X , 8F10« 1 ) 

READ (5.5000) NF 

5000 FORMAT (13) 

WRITE (6,1011) NF 

1011 FORMAT (lH0v34X.2I3) 

JUMP=NF 

IF (NF •EQ. 0) GO TO 5002 
R EAD ( 5.5) ( SPL ( I ) . I = 1 . NF ) 

READ (5,5)($LJE^( I),I==i,NF) 

WRITE (6.500) ( SPUI ) . I^1.NF) 

WRITE (6.500) (SLJET( I),I*1,NF) 

READ ( 5,5001) (COND( I ),I=1,12) 

WRITE (6.1012) 

1012 FORMAT (/) 

WRITE (6,1005) (COMDd ),I=1,12) 

5001 FORMAT (12A6) 

5002 continue 

IF (NF .EQ. 0) NF^NFO 
READ (5.5003) OASPL 

5003 FORMAT (FIO.O) 

WRITE (6,5004) OASPL 

5004 FORMAT ( 35X. FIO. 1 //) 

IF (OASPL .LE. 0.01) GO TO 501 
READ (5.5003) PS 
WRITE (6.5004) PS 
READ (5,3) (THACI ), I^lrNLlN) 

READ (5.3) (DAU ) .I=1.NL!N) 

WRITE (6,5005) (THA( I ) , 1=1 ,NLI N) 

5005 FORMAT (35X,5F10. 3) 

WRITE (6.5006) (D A( I ) . 1= 1 .NLl N ) 

5006 FORMAT ( IHO , 34X , 5 F10.3 ) 

END OF INPUT DATA READING 

501 CONTINUE 
DET=0. 

DL0H=0. 

00 1002 I=1.NLIN 
DO 1001 IFQ=1.NF 

1001 DDBN( I , IFQ)=0. 

DIFFE( I)=0. 

DIF»=L( I )=0. 

ET(1 )=H( I )«FP( I )/(C«12.) 

HG(I)=H(I) 

FPG( I)=FP( 1 ) 

BLOHGC I )=BLOH( I ) 

BL(I)=H(I)*BLOH(I ) 

DET=.DET+0.003’J^ET( I) 
DLOH=OLOH+0.0O3*BLOH( I ) 

1002 continue 

ITERAT=0 

ITEST=0 

10 INCR=0 

IF (NFSWCH. EQ. 0) GO TO 39 
GO TO 29 

11 CONTINUE 

IF (NFSWCH. NE. 0) GO TO 12 
DO 110 J=1.NLIN 
FP(J)=FPG(J) 
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110 5T( J )=H( J)^FP( J }/UZ^^C} 

12 CONTINUE 

DO 122 J=1,NLIN 
JR^U 

CALL ATTN( AM.EK J I.ODBOL) 

ODBPK=ODBOL*BLOH{ J) 

IF IET(J> .LT* ETAOO •AND. ET(J) .GT. 0.01) TP=PART 

IF (TR .EQ. 1. ) GO TO 115 

CALL OFFOPT {ET( J > ,TR,DOBPK,DDBPKO) 

ODBPK=DDBPKO 
115 CONTINUE 

PERUJ)^DDBOL 
PEAK! J )=:DDBPK 

C ALL BNDW ( DDBPK, ETC J ) ,FP U) , NF ,TR , DOB) 

NOTE THAT ODB IS USED HERE FOR TEMPORARY STORAGE. IT IS ALSO USED 
BELOW FOR PRINTOUT. 

00 121 IFQ=1.NF 

121 DDBNU, IFQ)=DDB( I FQ) 

122 continue 

IF (I TEST .GT. 0) GO TO 60 

ITEST IS A DUMMY VAPIASLE. WHEN IT IS .GT. 0 THE CALCULATION IS 
COMPLETE AND FINAL PRINTOUT IS PERFORMED BY THE GO TO STATEMENT. 
ITEST IS SET .GT. 0 BELOW. 

SUMOSQ^O 

DO 21 IFQ^ISTART, ISTGP 

totl=o. 

DO 20 J=1,NLIN 

20 TOTL=TQTL+OOBN( J, IFQ) 

DEV^GOALC IFQ)-TOTL 
DEVSQ = OEV=«'DEV»GOALUFQ) 

DEVSQ = OEV’«'OEV 

21 SUMOSQ=SUMDSQ+OEVSQ 
IF (LNGT--1) 31,41,41 


♦ NOTE THAT WHEN ANY ONE OF THE LOOPS BELOW 
*** IS FIRST REACHED THE NEEDED CALCULATIONS FOR 
THE DIFFERENCE SCHEME HAVE ALREADY BEEN 
MADE AT THE CENTRAL POINT. 


ETA LOOP 

29 LNGT=0 
1 = 0 

30 1 = 1 + 1 
1NCR=0 
GO TO 11 

31 INCR=INCR+1 

IF (INCR .EQ. 1) ETSAV=ET(I) 

GO TO (310,311,34), INCR 
31 e G=SUMOSQ 

ET( I)=ETSAV-'DET 
GO TO 312 

311 GMi=SUMOSQ 

ET ( n=ETSAV+DET 

312 FP ( I )r:ETC I )*C^12. /Hi I ) 

GO TO 11 

34 GP1=SUMDSQ 

CALL ANEWT ( ET SAV .DET .GMl .G.GPl.ETNEW > 
IF (ETNEW .LTo 0.) ETNEW=ABS ( E TNEW ) 
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39 


40 


41 


410 


411 


412 


540 


55 

C 

C 


C 

C 

C 

C 


56 

58 

60 

61 


65 


DIFFE( I)=ABS(ETSAV-ETNEW) 

ET<I)==£‘»'NEW 

FP( I )=ETNEW^C*12 . /Hi I } 

INCR=0 

IF (I •LI# NLIN) GO TO 30 
LOH LOOP 


LNGT=1 

1=0 


1 = 1 + 1 
INCR=0 
GO TO 11 
INCR=INCR+1 

IF (INCR .EQ* 1) SAVLDH=BLOH(l ) 

GO TO ( 410.411 ,41 2) » INCR 
G=SUMOSQ 

BLOH( I ) = SAVtOH-DLOH 

GO TO 11 

GM1=SUMDSQ 

BLOHC n = SAVLOH+DLOH 

GO TO 11 

GP1=SUMDSQ 

CALL ANEWT C SAVLOH,DLOH, GM 1, G , GPl . BLQHN I 
IF (BLOHN .LTo 0.) 8L0HN=ABS{ 8L0HN ) 
DIFFU I )=ABS(SAVLQH-BLOHN) 

IF LSAVLOH .LE* O.l) D!FFL{I)=0* 

BLOH< I )=BLOHN 
BL( I)=BLOHN*H( I ) 


INCR=0 


IF (I .LT. NLIN) GO TO 40 

CONTINUE 

ITERAT=ITERAT+1 

IF (ITERAT .EQ. 10) GO TO 61 


IF (ITERAT 
IF (ITERAT 
IF (ITERAT 
IF (ITERAT 
IF (ITERAT 

continue 


•EQ. 20) GO 61 

♦EQ. 30) GO TO 61 

♦EQ* 50) GO TO 61 

•EQ. 100) GO TO 61 
•EQ. 200) GO TO 61 


FRAC=0.01 
FR AC =0.05 
FRAC=0.02 
DO 56 J=1,NLIN 
TESTE=FRAC*£T( J ) 

TE STL=FRAC=(' BLOHC J) 


IF VARIABLES ARE NOT SUFFICIENTLY CLOSE TO RESULT OF PREVIOUS ITERATION, 
REPEAT THE ITERATION. 


IF (DIFFE(J) .GT. TESTE) GO TO 10 

IF (DIFFL(J) .GT. TESTL) GO TO 10 

CONTINUE 

ITEST=1 

GO TO 12 

CONTINUE 

DO 65 IFQ=1,NF 

ODBC IFQ)=0 

DO 65 J=1.NLIN 

DOB( IFQ)=008( IFQ) +OOBN(J, IFQ) 
WRITE (6,999) 
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999 FORMAT <1H1) 

WRITE (6f71) 

71 FORMAT {55X,23HLINER DESIGN PARAMETERS//) 

WRITE C6.68) (XPLNI I ) .1=1. 12) 

68 FORMAT (30X,12A6//) 

WRITE (6,710) 

710 FORMAT (60X.5H INPUT//) 

WRITE (6,711) 

711 FORMAT (45X.40H LINER H( INCHES) FP L/H//) 

DO 712 J=1,NLIN 

712 WRITE (6,713) J,HG( J),FPG( J) , BLOHG( J) 

713 FORMAT (49X.12.4X .F10.2.F10.0.F10.1) 

WRITE (6,998) 

998 FORMAT (IHO) 

WRITE (6.74) 

74 FORMAT (60X,18HCALCULATED RESULTS//) 

8LT=0. 

00 740 J=1.NLIN 
BL( J) =H( J)*BLOH(J ) 

740 BLT=BLT*BL( J) 

WRITE (6,741) 

741 FORMAT (38X, 5HL INER, 8X,3HL/H, 7X,2HFP ,4X,9HL( INCHES ) ,2X ,9HH ( INCHES) 
1,4X,3HETA//) 

DO 742 J=1,NL1N 

742 WRITE (6,743) J ,B LOH( J ) , FP ( J ) , 8L( J) ,H( J ) ,ET( J) 

743 FORMAT ( 40X . I 2 . 3X . F 1 0 . 2 . F 10. 0 . F 10. 1 .F 1 0. 2 . F 1 0. 3 ) 

WRITE (6,998) 

WRITE (6,744) BLT 

744 FORMAT (50X,24HT0TAL LENGTH (INCHES) = .FIO. 1) 

WRITE (6,998) 

NWHICH=0 

DC 7430 J=1.NL1N 

J=J 

IF (ET(J) .LE. ETAOO) NWHICH=J 
IF (NWHICH .NE. J) GO TO 7430 
WRITE (6,7431) NWHICH 

7430 CONTINUE 

7431 FORMAT ( 54X,5HLINER, I2,15H IS GFF-OPTIMUM) 

IF (NWHICH .EG. 0) GO TO 7439 

WRITE (6.7432) PART 

7432 FORMAT (54X,9HTHETA IS ,F4.1,20H TIMES OPTIMUM THETA//) 

7439 CONTINUE 

WRITE (6.75) 

75 FORMAT ( 55X,19HSP ECTRAL COMPARISON//) 

WRITE (6.76) 

76 FORMAT (56X,24HFREQ GOAL DOB//) 

DO 77 IFQ = I START, ISTOP 

77 WRITE (6.78) F ( IFQ) .GOAL ( I FQ) .00B( IFQ) 

78 FORMAT ( 50X,F10.0 ,2F 10.1 ) 

WRITE (6,79) 

79 FORMAT (IHO///) 

WRITE (6,80) I TER AT 

80 FORMAT ( 50X ,20HNUMBER ITERATIONS = ,14) 

IF (ITEST .NE. 1) GO TO 55 

CALCULATE PNL FOR HARO AND SOFT CONFIGUATIONS 

IF (JUMP .EQ. 0) GO TO 83 
CALL PNDBER ( 1 . NF .SPL . PNDBl, PNDB2) 

PNLHRD=PNDB1 
DO 81 IFQ=1,NF 
SPL( IFQ)=SPL( IFQ)-0DB( IFQ) 

IF (SPL(IFQ) .LT. SLJET(IFQ)) SPL( IFQ)=SLJET(IFQ) 
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81 CONTINUE 

CALL PN08ER ( 1 * NF .SPL . PN081. PNDB2) 

PNLSFT-PMOBl 
WRITE (6,999) 

WRITE C 6* 9999) 

9999 FORMAT ( IHO//// // //// ) 

WRITE (6,810) 

810 FORMAT ( 49X. 32HES TIMATEO PERCEIVED NOISE LEVELS//) 

WRITE (6,68) (XPLNd ) ,1 = 1, 12) 

WRITE (6.68) ( CON 01 1 ). 1=1. 12) 

WRITE (6.82) PNLHRD,PNLSFT 

82 FORMAT (51X,16HHARD WALL PNL - ,F10.1/ 

151X.16HS0FT WALL PNL = .FlO.l// 

248X,33HJET NOISE FLOOR EFFECT CONSIDERED) 

83 CONTINUE 

NOW SPECIFY BACKING DEPTHS, OPEN AREAS, ETC. 

IF (OASPL .LE. 0.01) GO TO 95 
DO 90 J=1,NLIN 

CALL CRITi2(ET( J) .AM.F IMP .RO . TG.E IG) 

FIMP=FIMP^ET( J) 

THETA( J)=REAL(FIMP) 

IF (ET<J) .LE. ETAOO) THETA( J)=PART«THETA( J) 

CHI(J)=AIMAG (FIMP) 

TH=THA(J) 

D=DAC J) 

CALL AR8D ( THET A( J) , CH I ( J ) , FP { J ) , S IGMA ( J ) , B( J ) ) 

90 CONTINUE 
WRITE (6.999) 

WRI^E (6,900) 

900 FORMAT ( 46X , 37HSP EC I F ICAT IONS OF LINER WALL GEOMETRY///) 

WRITE (6.68) (XPLNd >.1=1.12) 

WRITE (6,91) 

91 FORMAT (38X.5HLINER.8X.2HFP.7X.5HTHETA.6X.3HCHI.3X.i0HSIGMA(PCT> .2 
1X.6HB( IN.)//) 

DO 92 J=1,NLIN 

92 WRITE (6.93) J . FP ( J ) .THETA( J ) .CHI ( J) . S IGMA( J ) . B< J ) 

93 FORMAT ( 39X , 12 , 4X ,F10. 0, 2F10. 3 ,F10. 1 ,F 10.3 ) 

WRITE (6,998) 

WRITE (6.94) 

94 FORMAT (57X,14HF0R INPUT DATA//) 

WRITE (6,940) OASPL, PS 

940 FORMAT ( 52X .16H0A SPL IN-DUCT = .FiO.l/ 

152X, 19HPS(PSIA) IN-DUCT = , FlO.l//) 

WRITE (6,941) (THA( I ), 1=1.NLINER) 

941 FORMAT (45X .9HTH( IN) = .5F10.3//) 

WRITE (6,942) ( DA ( I ) , 1=1 , NL INE R ) 

942 format (45X.8HD(IN) = .5F10.3) 

95 CONTINUE 
N0CHMG=1 

IF (NOCHNG .EQ. 1) GO TO 1000 

STOP 

END 
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SUBROUTINE ANEWT (X,OX,GMl,G,GPl,XNEW) 

CALCULATE NEW GUESSED VALUE FOR ROOT, I.E, X, FOR AN EQUATION 
Y = FIX), FIND X SUCH THAT Y IS ZERO. 

IN PRESENT CASE F(X) IS THE FIRST DERIVATIVE OF G WRT X. 
APPLICATION IS FOR FINDING OPTIMUM (MAX OR MIN). 

METHOD USED IS NEWTON ITERATION. 

X(N+l)=X(N)-F(X(N)) /(DF(X(N) )/DX) 


THE DERIVATIVES ARE FOUND BY CENTRAL DIFFERENCE SCHEME. 

THUS THE CALLING PROGRAM MUST PROVIDE FUNCTIONAL VALUES AT THE POINT 
WHERE the DERIVATIVE IS DESIRED. AS WELL AS AT A STEP ON EITHER SIDE 
OF THE POINT. 

THESE VALUES ARE STORED IN GMl. G. GPl. 

X IS THE CENTRAL POINT. 

DX IS THE STEP SIZE. 

XNEW IS the CALCULATED NEW VALUE IN THE ITERATION FOR THE ROOT. 

DGDX=( GP1-GM1)/(2.*DX) 

DGDX2=XGP1-2.*G+GMI)/(DX*0X) 

F=DGDX 

0FDX=DGDX2 

XNEW=X-F/DFOX 

RETURN 

END 


SUBROUTINE ATTN { AM,ETA,DDBOL > 

CALCULATE PEAK ATTENUATION FOR A LINER, GIVEN ETA AND MACH NUMBER. 
NOTE THAT THIS IS VERSION 3 OF THIS PROGRAM 
MACH NUMBER EFFECTS ARE HEREIN HANDLED BY I NTERPOLAT I NG 
BETWEEN THE CURVES AT FIXED MACH NUMBERS. 


ARGUMENTS INPUT 

AM = FLOW-BY MACH NUMBER 
ETA = liner CHARACTERISTIC PARAMETER 
AGRUMENT OUT 

DDBOL = PEAK ATTENUAT ION/ ( L/H) , FOR OPTIMIZED RESISTANCE, 

LEAST ATTENUATED MODE RECTANGULAR DUCT. 

DIMENSION D(6,4),C(4),B{5) 

DIMENSION AMSTOR( 10) ,WQRK( 10) , ARG< 10) ,VAL( 10) 

DIMENSION ODBC 10) 

DATA { 8 ( K ) , K= 1 , 5 ) / . 41792 ,.57188,. 64086 , .239 05,14. 167/ 

DATA ( ( D ( I M , J ) , J= 1 , 4 ) , I M= 1 , 6 ) / 1 . 5475 , 4 5 . 4 97 , 1 8. 8 1 1 , 2 . 40 12 . 1 . 1 6 6 1 . 
110.937.9. 3986, 1.1 438.0.. 0. ,0. .0. .^2. 0987 .30. 741.12.255,2.4699,-3.5 
2123, 49.095,13.861 ,3. 52 51 ,-4.1 881 ,48. 918, 14.1 21 ,3.4267/ 

DATA (AMSTORI IM), 1M=1,6) /-.4.-.2,0. ..25..4, .5/ 

CALCULATE MACH=0 PEAK ATTENUATION. 
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EXP 1= ( ALOG ( ETA ) -A LOG ( B ( 2 D ) /8 ( 3 ) 

EXP2=-0*5»EXP1^EXP1 
TRM1=B<1)*EXP( EXP 2) 

TPM2= EXP (TP. Ml) 

COEFF=8(5)/(B(4)+ETA ) 

DDBNRM=C0EFF^TRM2 

IF (ETA oGE. 2.) GO TO 50 

IDENT=:0 

!TEST=0 

DO 25 IM=lf6 

IF (AM •EQ# AMSTOP(IM)) GO TO 23 
GO TO 25 
23 ITEST^l 

ident=im 

25 CONTINUE 

SET UP CURVES OF DELTA DB PK AT MESH POINT MACH NUMBEPS. 

LOOP ON lOUM IS FOR MACH NUMBER* 

THE J INDEX IS FOP THE PARAMETEPS IN THE FIT OF DELTA VERSUS ETA, 
AT THE selected MACH NUMBER. 

DO 31 I0UM=1.6 

IF (I test .EQ. 1) GO TO 26 

IM=IOUM 

GO TG 2? 

26 IM=IDENT 

27 continue 

IF ( IM .EO. 3) GO TO 29 

EXPN = D( IM,2)/(0( IM,4)+ETA)'-0{ IM,3) 

D E L = D ( I M . 1 ) * E ta**EXPN 
GO TO 30 

29 0£L=0. 

30 ODBC IM)=DDBNRM+DEL 

IF ( I test •EO. 1) GO TO 51 

31 continue 

interpolate on MACH NUMBER TO THE DESIRED VALUE. 

N0IM=:3 

ICOL=l 

IR0W=6 

C ALL ATSG (AM, AMS TOR . DDB ♦ WORK . I ROW • I COL . ARG • VAL * NO I M ) 

CALL LGRTRP (ARG, VAL,NDIM,AM,DDBOL) 

GO TO 52 

50 0EL=0. 

51 D 0 B 0 L = D D B N R M + D E L 

52 continue 
0DB0L=0D80L*0.4 

return 

END 


67 



o o r> n n r> o n o n n 


SUBROUTINE PN08ER ( I It 1 2t 06, PNOBlt PNDB2) 


THIS SUBROUTINE COMPUTES THE PNL FOR SPECTRUM 08, WHICH IS AN ARRAY WITH 
THE INDEX I CORRESPONDING TO THE NUMBERS OF THE THIRD OCTAVE BANDS* 

ARGUMENTS 

11 = FIRST BAND TO BE CONSICERED IN PNDB CALCULATION NORMALiy 1 

12 == LAST BAND 

OB ^ ARRAY CONTAINING THE SPECTRUM WHOSE PNL IS TO BE CALCULATED 
PN081 = PNL CALCULATED BY SMOOTHED CURVE FIT OF THE TABULAR 
DATA 

PNDB2 = PNL CALCULATED USING TABULAR DATA* 

DIMENSION DBK9), CB10(9) ,0B64( 15> ,DB50{1C) ,DB{ 27) 

DATA ( 061 CD, I = 1, 9) /63* 6 , 5 9* 7 ,56* 2 , 53 • ,50* 6,48., 45.6,44. 5, 43. / 

DATA C DBiO ( 1 ) * 1 = 1 .9 ) / 87. . 85. . 83. 1 .80. . 79.2 .78. .76. . 75. .74. / 
data (0664C I ), I =1,15)/ 112. ,111., 109. ,107., 106. ,105., 103. , 102. ,1 01. 

I, 100. , 1 ?0 . , 10 3 . , 1 00 . , 1 00 . , 98. / 

D A't A ( 0 B5 0 ( I ) . I =1 . 10 ) / 91 . . 89 . . 87 . . 86. . 8 6 . . 87 . . 88 . . 9 1 . . 94 . . 9 7. / 

SUMN 1=0.0 
SUMN2=0.0 
FMAX1=0.0 
MAX 2= 0.0 
DO 50 1=11. 12 
IF (I .GT. 15) GO TO 30 
£X = 0.0 30103’^{DB(I )-0B64CI ) ) 

PN1=64.0*10.0^*EX 

FN2=FN1 

IF (FNl .GT. 15. .OR. I .GT« 9) GO TO 40 
£X = <DB( I )-Obl( I )) /(DBIOC I )-DBl(I) ) 

FNLOW=10.0*^EX 

IF (FNLCW .LT. FNl) FN1=FNL0W 
FN2=FN1 
GO TO 40 
30 J=I-15 

DBC=D850( J ) 

EX = 0.030103*(DB( I )^DBC> 

FNI=50.0’«‘10.0^»EX 

IF (D8(I) .GE. OBC) GO TQ 35 

X=( OBC-DBC 1 ) 

FNI=FN1’«'<1. 0+0.02 79»X/( 1.0+0. 6«X) ) 

FN2=FN1 
GO TO 43 

35 X = DB(n-DBC 

ARG=D.2^3.1415927«X 

FS=SIN(APG) 

COR=FS 

IF CX .GT. 10.0) GO To 37 
CnR=< 3.5-0.275*X) ^FS 

IF (X .LT. 3.0) COR=COR+O.6=«^X’^(3.0-X) 

37 FN2=FND( 1 .0+0. 01 ♦COR ) 

40 $UMN1=SUMN1+FN1 
SUMN2=SUMN2+FN2 
IF (FNl .GT. FMAXl) FMAX1=FN1 
IF, <FN2 .GT. FMAX2) FMAX2=FN2 
50 continue 

F 1 =0 . 8 5«F MAX 1+0 . 1 5^ S UMNl 
?"2=0.85*FMAX2+0.15*SUMN2 
PM OB 1 = 10 . 0* ( 4. 0+A LOG 1Q(F1)/0.3C103) 

PNDB 2=1 0 . 0 ♦ ( 4 o 0 + A L DG 1 0 ( F 2 ) / 0 . 3 0 1 0 3 ) 

RE'^URN 

END 


68 



o n o o r> n n o o o n o o o o n o o r> o n 


SUBROUTINE BNOW ( DDBPK , ETA, FP , NF , TR, DOB) 


CALCULATE ATTENUATION SPECTRUK FOR LINER WITH SPECIFIED PEAK 
ATTENUATION, ETA, FREQUENCY OF PEAK, AND RESISTANCE RATIO. 

ARGUMENTS IN ' 

DDBPK - LINER PEAK ATTENUATION, MUST BE ADJUSTED EXTERNALLY IF 
LINER IS OFF-OPTiMUM. 

ETA - LINER CHARACTERISTIC PARAMETER, AT PEAK FREQUENCY. 

FP - LINER PEAK FREQUENCY 

NF - NUMBER OF THIRD OCTAVE FREQUENCIES FOR CALCULATING SPECTRUM 
TR - RATIO OF RESISTANCE/OPTIMAL RESISTANCE. 

ARGUMENT OUT 

9DB I) - ARRAY CONTAINING THIRD OCTAVE ATTENUATION SPECTRUM 
COMMON /FREQ/ F(27) 

DIMENSION ET(5) ,THR 1 { 5 ,5 ) , A1 ( 5 ,4) , BH 5 ,4 ) , Yl (5) 

DIMENSION WORK! 5) ,ARG( 5) ,VAL( 5) 

DIMENSION A20PT(3),B20PT(3),THR2<4),A2(3),B2<3),BPWL(2),0DB<27) 

DATA (ET( I ),I=1,5) / .1 , . 4 , 1. , 2. , 10. / 

DATA ((THBld.J), J-1.5).I=1,5) /. 1 ,. 69,1 . , 3. 05 , 10. , . 1 , .6* 1. . 3. 7, 1 
10., . 1 , . 55 , 1. , 2. 8 , 10. , . 1 , . 57 , 1. , 2 . 42 , 10. , . 1, .484 , 1 . , 2. 25 , 10. / 

DATA ( (Al( I,J) ,J=1,4) ,1=1,5) /-. 051 036,-. 036, . 0 84, -.119, -.086, 
1-. 086, .262, -.276,-. 215,-. 215, .105, -.416,-. 314, -.314, -.037, -.553 
2, -.444, -.444, -.143/ 

DATA ((BK I,J),J=1,4),I=1,5) /Q. ,0.089, -.103, -.352 , 0. , . 149 ,-. 348, 
1-.961. 0.,. 235,-. 417,-1. 132, 0. , .415.-.471 ,- 1. 192 , 0. , . 346 ,-.4C 1 ,- 
21.255/ 

DATA <A20PT(I), 1=1,2) /.316,.328/ 
data (B20PT<I).I=1,2) /.083,.025/ 

DATA (THR2( J), J=1 ,4) /.I , 1. , 1. 9, 10. / 
data (A2(J),J=1,3) /.312.. 312,.236/ 

DATA (B2( J ) ,J=1,3 ) /O. ,.833.1.107/ 

DATA NETA,NSEG1,NSEG2 /5,4,3/ 


CHECK TO SEE IF ETA EQUALS ONE OF THE STANDARD VALUES, FOR FI, 

1WHICH=0 
T0L=.005 
DO 10 1=1, NET A 
I = I 

ETLD=ET( I)*( l.-TOL) 
eTHI=ET( I )«( l.+TOL) 

IF (ETA ,GE. ETLO .AND. ETA .LE. ETHI ) IWHICH=I 
IF (ETA .GT. ETCNETA)) IWHICHsNETA 
10 CONTINUE 

CALCULATE QUANTITIES FOR FI AT EACH STANDARD ETA WHEN INPUT ETA IS 
NOT A standard VALUE. 

DO 40 I=1,NETA 
JWHICH=0 
DO 20 J=1,NSEG1 
J = J 

IF (TR .gE. THRKI.J) .and. TR .LE. THRl ( I . J+1 ) ) JWHICH=J 

20 continue 

IF (JWHICH .GT. 0) GO TO 30 
IF (TR .LT. THRKI.D) JWHICH=1 
IF (TR .GT. THRK I,NSE61) ) JWHICH=NSEG1 
30 Yl( I )=A1(I .JWHICH )+Bl( I. JWHICH)*AL0G10(TR) 

40 CONTINUE 
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IF Cl WHICH •GT. 0) GO TO 45 
INTERPOLATE THE VALUE FOR NON-STANDARO ETA. 

NDIM=3 

CALL ATSG (ETA, ET, Y1 , WORK, NETA,1 ,ARG,VAL,NDIM) 

CALL LGRTRP (ARG,VAL,NDIM,ETA,Y1ETA) 

GO TO 46 

45 Y1ETA-Y1( IWHICH) 

46 CONTINUE 
F1QFP=10.«»Y1ETA 

MOW calculate F2/FP, DO OPTIMAL CASE FIRST. 

IF (ETA .GTa 1.581 GO TO 60 

JWHICH^l 

GO '*’0 61 

60 JWHICH-2 

61 CONTINUE 

Y20PT=A20PT(JWHICH)+B2QPT( JWHICHI^ALOGIOCETA) 

FR20PT = 10.^>«'Y20PT 
IF (Tr oME. 1.) GO TO 65 
F20FP=FR20PT 
GO TO 85 
65 CONTINUE 
JWHICH=0 
00 70 J=lfNSEG2 
J=J 

IF (TR .GE. THR2(J) .AND. TR .LE. THR2(J+1)> JWHICH=J 
70 continue 

IF (JWHICH .GT. 3> GO TO 80 
IF (TR .LT. THR2(D) JWHICH=1 
IF (TR .GT. THR2(NSEG2>) JWHICH=NSEG2 
80 Y21=A2( JWHICH) +82 (JWHICH)«ALQG10(TR) 

F0FF2=10.^*Y21 
F 2 0 FP = FQ F F 2^^^ ( F R 20 PT / 2 • 05 ) 

85 CONTINUE 

F20FP=F20FP^1.34/ 2.15 

CALCULATE OFF-PEAK DELTA PWLS. 

8 P WL ( I ) = A LOG 10 ( . 5 ) /A LOG 1 0 ( F 10 F P ) 

BPWL( 2)=AL0G10( .5) /ALOGIO ( F20F P) 

00 90 I=1,NF 
J = 1 

FR=F( I ) /FP 

IF ( FR .GT • 1. ) J *2 

Y = BPWLC J)’^'AL0G10( FR ) 

FRAC=10.«*Y 
90 DDB( I )==DOBPK*FRAC 
RETURN 
END 
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SUBROUTINE OFFOPT {ETA.TR.DDBPK.DOBPKO) 

CALCULATE PEAK ATTENUATION OF A LINER FOR WHICH THE RATIO OF 
RESISTANCE/OPTIMAL RESISTANCE IS NOT 1. 

ARGUMENTS IN 

ETA - VALUE OF FH/C AT FREQUENCY OF PEAK ATTENUATION 
TR - RESISTANCE/OPTIMAL RESISTANCE, AT FREQUENCY OF PEAK ATTENUATION 
DDBPK - PEAK ATTENUATION, OPTIMAL 
ARGUMENT OUT 

DDBPKO ^ PEAK ATTENUATION, OFF OPTIMAL 

DIMENSION ET(3) .AC3.3) »8C3*3) ,THRC3.4) •YG(3) 

DIMENSION WORKC5) ,ARGC5),VALC5> 

DATA (ETU), 1=1*3) /.I*. 4.1*/ 

DATA ( { A ( I * J ) ♦ J= 1 ♦ 3) * 1 = 1 , 3 ) /O* , 0* ,**• 151 , Q* ,0* j*”* 116, 0® ,0*f"*»110 
1 / 

DATA <<BU.J)*J=1*3).1 = 1*3) /- 1. . * 382 • . 753 . “1 . * . 790 . 1. 167 . -1...8 
175,1.248/ 

DATA ( (THR(I,J ) ,J=1,4) , 1=1,3) /. 1 , 1. , 2. 55 , 10« , • 1 ,1 . , 2. 03 , 10^ , •! 

1* 1 . * 1 • 97* 10* / 

DATA NETA,NSEG /3,3/ 

IF RESISTANCE RATIO IS EQUAL TO OR LESS THAN 1. USE SINGLE LINE 
FOR ALL ETA. 

IF <TR .GT* 1.) GO TO 20 
y=A ( 1, D+B a, D^ALOGIOITR) 

GO TO 60 

CALCUALTE BETAS AT STANOARO ETAS 

20 IWHICH=0 

DO 30 I=1,NETA 
I = I 

IF (ETA *£0* ETm> IWHICH=I 
30 CONTINUE 

IF (ETA .GE. ET(NETA)) IWHICH=NETA 

DO 50 1=1, NET A 

JWHICH=0 

DO 40 J=2*NSEG 

J = J 

IF (TR *GE. THR(I.J) .AND. TR *LE. THR(I*J+1)) JWHICH=J 
40 CONTINUE 

IF (TR *GT. THR(I,NSEG)) JWHICH=NSEG 
50 YG(1 ) =An* J WHICH) +B( I . JWHICH) *AL0G10(TR) 

INTERPOLATE THE VALUE FOR NON-STANDARD ETA. 

IF (IWHICH *EQ* 0) GO TO 55 
Y=YG( IWHICH) 

GO TO 60 
55 CONTINUE 
ND1M=^3 

CALL ATSG (ETA, ET.YG, WORK, NETA, 1 ,ARG,VAL,NDIM) 

CALL LGRTRP (ARG, VAL,NDIM,ETA,Y) 

BETA0=10.**Y 

DDBPKO=DDBPK/BETAO 

RETURN 

END 
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C 

C CALCULATES BRANCH PCINT IMPEDANCE R* THETA. FIRST.SECONO MODE 
C ^ 

S U BROUT INE CR I T 12 ( ET A , FMA C H, F I PP , R , T , E I G ) 

COMPLEX FIMP.EIG 
IF(FMACH#LT#0.0) GO TO 20 

RAT = I.O-0.276^FMACH*(1«O+O.8^FMACH)/CETA’8'*2+O#8’«'FMACH> 
IFCETA.LT.1.0) GO TO 10 
DT - 3.2^FMACH^(3.5-ETA) 

IF(DT.LT»0*0) DT=0*0 
GO TO 40 

10 DT = 30.0^FMACH^CETA«C 2#0-ETA>-0.7> 

GO TO 40 

20 RAT * 1.0-FMACH*l0#13-(2*75+2.5*FMACHI4ALOG10<ETA) ) 
IF(RAT.LT.1.0) RAT=i.O 
IFCETA.LT.1.0) GO TO 30 
OT = 2.8=5tFMACH’e^(7#0-4.04ETA) 

1F(OT.GT.O.O> OT^^O.O 
GO TO 40 

30 DT = 40.0*FMACH4:(ETA^( 1.6-ETA l-0#39) 

40 R = RAT^2.39 
T * 28.0+DT 
ANG = T/57.29578 
FI ^ R^COSIANGI 
F2 - R^SINCANG) 

EIG = CMPLXCF1,F2) 

IASYM=0 

CALL 1 RT 0 S 1 ( E I G , E TA , F MAC H i F I M P 1 1 A S YM ) 

RETURN 

END 


CALCULATES Z/ (RHO«C*ETA ) RECT. DUCT,D I SP.CONT, BCTH WALLS SAME IMPED. 
DUCT WALLS AT Y = -AO.'«-AO 

S U 8R 0 U T I N E ZR T D S 1 < E I G . £ T A » F M A C H . F I M P . I A S YM ) 

COMPLEX EIG.FIMP.UNCP.StQ.CPl 
UNCP == CMPLXIO.O, 1.0) 

PI = 3.1415927 

FMl = 1.0-FMACH*’«'2 

CPI ^ FM1*(EIG/PI/ETA)**2-1.0 

S = CSQRTCCPl) 

Q = ( 1.0+UNCP*FMACH*S)/FM1 

F IMP = UNCP’9‘PI^CC0S( EIG) ^Q«^2/ (EIG^CS INI EIG) ) 

IFUASYM.EQ.O) GO TO 10 

FIMP --UNCP«PI^CS1N(EIG)’«'Q^*2/CEIG*CC0S<EIG) ) 

10 CONTINUE 
RETURN 
END 
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SUBROUTINE ARBD rTHETA , CHI , FP , SIGMA, B) 


C 

C CALCULATE AREA RATIO ANO BACKING DEPTH FOR LINER DESIGN^ 

C IMPEDANCE MODEL USES ADDITIVE FLOW-BY CORRECTION* 

C UNLESS NBYSW IS OTHER THAN ZERO, E.G* 1 

C 

C ARGUMENTS 

C theta = NORMALIZED RE^SISTANCE 

C CHI = NORMALIZED REACTANCE 

C FP = LINER TUNE FREQUENCY 

C SIGMA = PERCENT OPEN AREA 

C B = CAVITY BACKING DEPTH, INCHES 

C IN COMMON 

C PS ^ DUCT STATIC PRESSURE* PSI 

C AM = DUCT MACH NUMBER 

C T = duct TEMPERATURE, DEGREES RAMKINE 

C OASPL = DUCT OA5PL 

C TH ^ FACING SHEET THICKNESS, INCHES 

C 0 = FACING SHEET HOLE DIAMETER, INCHES 

C NBYSW = SWITCH TO TELL DESIRED IMPEDANCE MODEL TO USE* 

C IF ZERO, USE ADDITIVE FLOW-BY 

C IF HON-^ZERG, USE LINEAR FLOW--8Y. 

C 

COMMON /PESANT/ P S, AM , T, OASPL , TH, D 
COMMON /FL08Y/ NBYSW 

NAM EL I ST /NAl / PS , AM , T ,OASPL , TH, D, THETA , CH I , FP , OMEGA 

NAMELIST /NA2/ RHO,C,Z 

NAMELIST /NA3/ AMU*GNU*SP 

AMSAV=AM 

AM=A8S( AM) 

0MEGA=2**3*141592 7^FP 

RHO=2.710982^PS/T 

C=49*02*SQRTIT) 

VINF=ABS( AM)^C 

AM ( 7 . 30248E- 7 ) « I T*^l • 5 > / < T+ 1 98* 72 ) 

GNU=AMU/RHO 

PREF=2.*14*5038«144*E-10 
SP=PREF«10*^=^ (GASPL/20# ) 

BETA=< l*+TH/0 )*saRT( 8**GNU=«'GMEGA) /C 
ZETA=0*3^ABS(AM) 

Z = SQRT ( THETA** 2+C HI><'*2) 

GAMMA=SP*32#17/(RH0*C*C*Z) 

GAMPRM = *424*SQRT< 2. )*GAMMA 
IF (NBYSW .NE* 0) GO TO 10 
FIRST=3ETA + ZETA 

SEC0ND=FIRST*FIRST+4**THETA*GAMPRM 
GO TO 11 

10 FIRST=BETA 

S EC0ND=BETA*B£ TA+ 4**THETA*GAMMA* ( 1. +0.006* VI NF) 

11 SIGMA={FIRST+SQRT(SEC0ND) )/(2**THETA) 

DELTA=0*85*D*( 1.-“.7*SQRT( SIGMA )) /( l.+305.*( A8S (AM) )**3 ) 
ALE=TH+OELTA 

ARG=0MEGA*ALE/( SIGMA*C*12. )~CHI 
ARG=1./ARG 

b=atan( arg)*c*12. /omega 

IF (B .LE. 0.) B=B+3.1415927*C*12. /OMEGA 

sigma^sigma*ioo. 

AM^AMSAV 

RETURN 

END 
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Printed Output 


The initial p^e of the output section is a listing of sample data in proper format 
for running the program. The remaining p^es contain the liner design information 
calculated by the program, Including a comparison between the goal attenuation spec- 
trum and the predicted spectrum for the liner design. The variable names used are, 
as nearly as possible, consistent with the symbols in the body of the report. Obvious 
exceptions to this are the Greek symbols, such as rj, which are not available on the 
computer printer are therefore spelled out, such as ’*ETA”. 


TABLE XK. - SET OF INPUT DATA - LINER DESIGN 

[Note that initial guesses must be made for tune frequency and liner lengths to begin iteration 
solutions. ] 

£•00 0.70 

1 2 

0*40 519.0 

500. 2500. 

2.0 4.0 

9.70 9.70 

8 22 


1.0 

4.7 

7.5 

20.1 

10.1 

10.0 

14.5 

9.3 

15.1 

15.0 

14.5 

13.2 

13.3 

8.9 

3.1 


69.2 

70.2 

70.6 

70. 7 

70.3 

69.9 

70.1 

71.3 

73.6 

75.9 

88.5 

78.5 

78.4 

82.2 

75.2 

78.2 

75.4 

73.0 

7 3. 7 

70.9 

67.4 

63.5 

61.4 

55.0 

69.2 

70.2 

70.6 

70.6 

70.1 

69.3 

68. 1 

66.9 

65.5 

64.0 

62.4 

60 .8 

59.1 

57.4 

55. 8 

53.9 

52.2 

50.4 

48.3 

46. 1 

43. 6 

40.6 

36.9 

32.5 

AT 

500. 

FE5T SIDE 


120. DEGREES 





144.0 

15.0 

1.500 0.020 

0.250 0.050 
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TABLE XX. - LINER DESIGN PARAMETERS 


[Liner for Hamilton standard based on noise estimate, exhaust, and 
90.0-PNdB goal, ] 




INPUT 





LINER 

H( INCHES) 

FP 

L/H 



1 

9.79 

500. 

2.0 



2 

9.70 

2500. 

4.0 




CALCULATED RESULTS 



LINER 

L/H 

FP 

LMNCHES) 

H( INCHES) 

ETA 

1 

2.35 

502. 

22.8 

9.70 

0.363 

2 

6.24 

3006. 

60.6 

9.70 

2.176 


TOTAL LENGTH (INCHESI = 83.4 


LINER 1 IS OFF-OPTIMUM 

THETA IS 2.0 TIMES OPTIMUM THETA 

SPECTRAL COMPARISON 


FREQ 

GOAL 

DOB 

250. 

1.0 

5.1 

315. 

4.7 

7.1 

400. 

7.5 

10.1 

500. 

20.1 

14.2 

630. 

10. 1 

12.9 

800. 

10.0 

12.0 

1000. 

14.5 

11.7 

12 50. 

9.3 

11.8 

1600. 

15.1 

12.6 

2000. 

15.G 

13.8 

2500. 

14.5 

15.6 

3150. 

13.2 

15.9 

4000. 

13.3 

9.7 

5000. 

8.9 

6.2 

6300. 

3.1 

3.9 


NUMBER ITERATIONS = 6 


TABLE XXI. - ESTIMATED HARD AND SOFT 
WALL PNL'S AT 152. 5- METER (500- FT) 

SIDELINE AND 120° 

^ [Liner for Hamilton standard based on noise 
estimate, exhaust, and 90.0-PNdB goal. ] 

HARO WALL PNL = 100.5 

SOFT WALL PNL * 89.5 

JET NOISE FLOOR EFFECT CONSIOEPED 



TABLE XXII. - SPECIFICATIONS OF LINER WALL GEOMETRIES 


INCLUDING TUNING, IMPEDANCES, INPUT WALL THICKNESSES, 
AND HOLE DIAMETERS 


[Liner for Hamilton standard based on noise estimate, exhaust, 
and 90.0~PNdB goal. ] 


LINER 

FP 

THETA 

CHI 

SIGMA(PCT) 

B( IN. ) 

1 

502. 

0.595 

0.191 

23.7 

2.779 

2 

3006. 

1.115 

-0.794 

12.0 

0.541 



FOR INPUT DATA 




GASPt 

IN- DUCT = 

1441. 




PS(PSIA) IN-DUCT 

= 

15.0 



THUN) = 1.500 0.020 

DUN) = 0.250 C.050 


Definition of variables: 

FP Peak frequency 

THETA Specific resistance 

CHI Specific reactance 

SIGMA Face sheet percent open area 

B Cavity backing depth 

TH Face sheet thickness 

D Face sheet perforation diameter 
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Figure 2. - Suppressor design procedure. 



Figure 3. - Composite noise prediction spectrum for fan at 30. 5-meter QOO-ft) sideline, 
120° from inlet direction. Aerodynamic parameters given in table I. 
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Specific resistance, 



Frequency, f. kHz 


Figure 13. - Adjusted attenuation estimate compared with measured spectra 
at several speeds from QF-IA fan inlet tests from reference 18. 



Specific reactance, x 


Figure 14. - Illustration of equal attenuation contours and liner impedance model for three 
different frequencies. (Circles are approximations to actual curves. ) 
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Frequency parameter, t?^ Hfp/ c 


Figure 17. - Bandwidth characteristics for optimum resistance liner. 
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(a) Back cavity depth. 



Frequency parameter. 7 } = Hfp/c 


(b) Open area ratio. 

Figure 20. - Chart for liner wail selection for optimum imped- 
ance. Duct Mach number, M, 0.4; duct temperature, T. 

288 K (519° R); static pressure, p, 1.013 N/m^ (14.7 psia); 
overall sound pressure level, OASPL, 140 dB; sheet thick- 
ness, t, O.OSlcentimeter (0.020 in, 1; hole diameter, d, 

0. 127 centimeter (0. 050 in. ). 


Attenuation, APWL, AdB 



Frequency, f, kHz 


Fig u re 21. - Meas u red a nd estimated acoustic line r performance. Gene ral 
Electric Quiet Engine C inlet liner; 80 percent speed. (NASA unpublished 
acoustic tests. ) 


NASA- Langley, 1976 E-8428 


91 






















